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INTRODUCTION 


The  use  of  fiber-reinforced  plastics  (FRP)  is  increasing  in  marine  structural  applications 
ranging  from  naval  ships  (e.g.,  hulls  of  minehunters,  deckhouses,  appendages)  to  high-performance 
yachts  to  offshore  structures  (e.g.,  piers,  risers,  topside  structures).  This  trend  is  expected  to 
continue  for  two  reasons.  First,  the  price  of  FRP  material  is  steadily  decreasing,  approaching  a 
level  comparable  to  that  of  metals.  Second,  low-cost  and  high-quality  manufacturing  techniques, 
such  as  autoclave  and  resin  transfer  molding  processes,  have  been  developed,  making  production 
of  FRP  more  attractive. 

FRP  applications  in  marine  environment  often  prompt  questions  about  fatigue  of  the 
material  because  ship  structures  are  constantly  subjected  to  cyclic  wave-induced  loading.  In 
addition  to  wave  action,  submersible  vehicles  can  experience  a  sufficient  number  of  diving  and  re¬ 
surfacing  cycles  to  cause  low-cycle  compressive  fatigue  damage  [Smith  (1990)].  Fatigue  issues  in 
FRP  marine  structures,  have  received  little  attention  because  the  design  of  most  FRP  ships  has 
been  constrained  by  the  material  stiffness  rather  than  by  strength  limitations.  Such  designs  have 
resulted  in  low  operational  stress  levels  at  which  fatigue  failures  were  not  expected  during  the 
specified  service  life.  This  situation  is  changing,  however,  due  to  two  main  factors.  First,  when 
some  FRP  ships  complete  the  operational  period  initially  intended,  a  question  may  arise  about 
whether  this  period  can  be  extended.  Second,  FRPs  are  increasingly  preferred  for  replacement 
weight-saving  solutions,  in  which  case,  many  critical  parts  can  experience  fatigue-stimulating  stress 
levels.  In  addition  to  these  two  factors,  there  is  always  the  issue  of  fatigue  in  the  structural  joints, 
which  are  usually  zones  of  stress  concentration,  complex  multi-axial  loading,  and  low-quality 
workmanship.  Therefore,  fatigue  is  becoming  an  important  issue  for  FRP  marine  structures. 
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Design  against  fatigue  failure  is  sometimes  difficult,  even  for  structures  made  of  metals 
such  as  steel  and  aluminum,  because  of  the  interrelated  influence  of  a  number  of  factors.  These 
include  microstructure  and  material  properties,  dimensions  and  geometry,  production  techniques, 
loading  history,  and  environmental  effects.  Nevertheless,  design  against  fatigue  in  metal  structures 
has  been  relatively  successful;  because,  regardless  of  the  effects  of  these  factors,  metals,  as  well  as 
welded  joints,  tend  to  follow  power  law  stress-fatigue  life  relations  and  the  linear  fatigue  damage 
accumulation  [Palmgren  (1924),  Miner  (1945)].  By  contrast,  design  against  fatigue  in  FRP 
structures  encounters  two  major  difficulties.  First,  the  number  of  factors  affecting  their  fatigue  life 
inevitably  increases  given  the  complex  structure  of  the  FRP  laminate,  which  is  a  combination  of 
two  components,  matrix  and  fibers,  with  quite  different  properties.  The  second  difficulty,  and 
probably  more  important,  is  that  the  limited  number  of  studies  for  marine  applications  suggest  that 
the  damage  accumulation  processes  in  FRP  may  not  obey  the  simple  linear  relations  encountered  in 
metals.  Such  a  difficulty  suggests  that  current  fatigue  guidelines  for  steel  ship  structures  may  be 
totally  inadequate  for  FRP  ship  structures.  A  simple  example  is  the  fact  that  FRP  damage 
accumulation  may  be  affected  by  the  loading  frequency,  which  is  not  usually  the  case  with  metals; 
this  means  that  a  frequency  parameter  must  be  included  in  the  damage  accumulation  model  for  an 
FRP  structure. 

In  order  to  address  the  aforementioned  difficulties,  an  experimental  and  analytical 
investigation  was  undertaken  at  the  Center  for  Structural  Djmamics  Research  of  The  George 
Washington  University.  The  investigation  concentrated  on  two  major  tasks.  The  first  task  was  an 
extensive  literature  search  for  publications  that  concern  fatigue-related  issues  for  FRP  composite 
laminates  and  joints.  Specific  attention  was  paid  to  all  pertinent  fatigue  damage  accumulation 
theories  proposed  in  the  literature.  The  second  task  encompassed  an  extensive  testing  program  to 
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enhance  the  understanding  of  fatigue  processes  in  FRP  laminates,  provide  enough  experimental 
data  for  assessment  and  improvement  of  the  methodologies  identified  under  the  first  task,  and  to 
serve  as  a  basis  for  development  of  new  such  methodologies. 

This  report  summarizes  the  most  important  findings  under  both  tasks.  It  contains  two 
major  sections,  each  dealing  with  one  of  the  tasks.  The  first  section  describes  the  physical  aspects 
of  fatigue  damage  in  FRP  laminates,  the  analytical  means  of  characterizing  constant  amplitude 
fatigue  data,  and  the  more  significant  damage  accumulation  studies  reported.  The  issues  related  to 
the  effect  of  parameters  such  as  loading  frequency,  void  and  fiber  content,  temperature,  and 
environment,  (seawater  and  humidity)  are  not  considered  here.  These  efforts  are  recommended  for 
future  studies.  The  second  section  discusses  the  testing  program  and  the  preliminary  results 
obtained. 


FATIGUE  DAMAGE  OF  FRP  LAMINATES 

PHENOMENOLOGICAL  ASPECTS 

In  FRP  laminates,  fatigue  damage  may  manifest  itself  as  either  fiber  breakage  or  matrix 
cracking,  but,  in  most  cases,  it  is  a  combination  of  both.  The  particular  form  of  damage  depends 
on  the  laminate  reinforcement  type  (uni-,  bi-  or  multi-directional),  the  loading  mode  (tensile  or 
compressive),  and  the  orientation  of  the  loading  with  respect  to  the  fibers. 

When  a  unidirectional  composite  is  loaded  by  longitudinal  tensile  stress  parallel  to  the 
fibers,  the  surrounding  matrix  may  undergo  a  fatigue  process  of  cracking  in  the  direction  transverse 
to  the  fiber  direction  before  any  local  fiber  breaks  [Talreya  (1987)].  These  cracks  are  initially 
randomly  distributed  and,  if  the  local  strains  are  low,  remain  arrested  by  the  surrounding  fibers. 
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When  the  local  strains  become  higher  than  a  certain  threshold,  a  fiber  may  break;  the  event  being 
assisted  by  the  previously  created  matrix  cracks.  A  broken  fiber  causes  shear  stress  concentration 
and  debonding  of  the  fiber-matrix  interface  near  the  two  broken  tips.  The  debonded  zone  acts  as  a 
stress  concentrator  leading  to  additional  matrix  cracking  in  a  direction  perpendicular  to  that  of  the 
loading  and  the  fiber.  Such  progressive  crack  growth  eventually  leads  to  a  visible  macro-crack 
extending  to  the  free  interface  of  the  laminate.  Unidirectional  composites  are  known  to  be  quite 
sensitive  to  loading  inclined  to  the  direction  of  the  fibers:  only  a  few  degrees  of  misalignment  can 
cause  the  strain  at  failure  to  reduce  by  half  [Reifsnider  (1990)].  In  such  a  case,  the  predominant 
damage  mode  is  matrix  cracking  along  the  fiber-matrix  interface;  this  is  the  process  that  controls 
the  fatigue  strength. 

Unidirectional  composites  are  rarely  used  in  practice.  More  common  are  the  angle-ply 
laminates  in  which  unidirectional  laminas  are  stacked  in  alternating  directions.  Angle-ply  laminates 
are  usually  composed  by  0/90°  stacking  (cross-ply  laminates),  although  plies  at  other  angles, 
usually  ±  45°,  are  also  included  in  the  stack  to  obtain  more  uniform  (isotropic)  material  properties 
in  all  directions.  The  primary  fatigue  damage  of  angle-ply  laminates  consists  of  matrix  cracking 
along  the  fiber  direction  in  those  plies  whose  fibers  are  not  aligned  with  the  loading  direction 
[Reifsnider  (1980)].  The  cracking  process  continues  until  each  off-axis  ply  attains  a  stationary 
crack  saturation  state.  After  that,  no  new  cracks  develop  because  adjacent  cracks  reduce  the  local 
stress  to  such  an  extent  that  the  material  between  them  cannot  develop  new  cracks.  Such  a 
saturation  state  is  characteristic  for  the  matrix  and  fiber  materials  and  the  laminate  stacking 
sequence  but  seems  to  be  independent  of  the  loading.  Therefore,  it  is  commonly  referred  to  as  the 
“characteristic  damage  state”  (CDS)  [Talreya  (1987)]. 

Fatigue  damage  is  also  assisted  by  interlaminar  shear  stresses  and  normal  stresses  in  the 
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thickness  direction.  Such  stresses  naturally  exist  around  free  edges  of  undamaged  laminates  but 
also  develop  in  areas  damaged  to  the  point  of  CDS.  Therefore,  the  damage  state  following  CDS  is 
characterized  by  ply  delamination.  Delamination  further  promotes  fatigue  damage  by  allowing 
cracks  in  adjacent  plies  to  connect  with  one  another  in  the  through-thickness  direction.  Once  the 
matrix  is  weakened  enough  by  cracking,  fibers  start  to  break  and  ultimate  failure  occurs  at  the 
maximum  load  of  a  cycle. 

In  addition  to  angle-ply  laminates,  two  other  types  of  FRP  have  found  a  wide  application  in 
marine  structures.  These  are  the  chopped  strand  mat  (CSM)  and  woven  roving  (WR)  laminates. 
CSM  composites,  which  are  common  in  commercial  marine  craft,  do  not  have  a  continuous, 
oriented  fiber  system:  the  fibers  are  relatively  short  and  are  randomly  dispersed  in  the  matrix  during 
manufacturing.  It  has  been  found  that  the  fatigue  damage  they  develop  is  similar  to  those  of 
laminates  with  continuous  fibers  [Owen  and  Howe  (1968)].  The  matrix-fiber  interface  first 
debonds  along  those  fibers  that  happen  to  be  oriented  perpendicular  to  the  loading  direction.  As 
the  load  cycling  continues,  damage  spreads  to  fibers  aligned  at  smaller  angles  to  the  loading 
direction.  Damage  accumulation  seems  to  be  accelerated  by  matrix  cracking  in  the  zones  where 
the  fiber  density  is  lower  as  a  result  of  the  manufacturing  process.  This  type  of  damage  usually 
governs  the  fatigue  behavior  of  CSM  laminates  at  low  stress  cycles,  which  suggests  that  fatigue 
properties  of  these  laminates  in  marine  environment  may  be  predominately  matrix-controlled. 

WR^  laminates  are  common  on  larger  naval  ships  such  as  minehunters.  Three  stages  of 
fatigue  damage  are  usually  observed  when  a  WR  laminate  is  loaded  in  the  warp  direction 


‘  A  woven  fabric  consists  of  two  sets  of  interlaced  yams  in  two  (usually  perpendicular) 
directions;  one  of  the  directions  (usually  in  the  direction  of  the  loading)  is  referred  to  as  ‘warp’  and 
the  other  is  termed  as  ‘weft’.  The  crossovers  are  known  as  ‘undulations’. 
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[Tanimoto  and  Amijima  (1975)].  First,  transverse  cracks  initiate  in  the  matrix  at  the  undulations 
and  gradually  reach  a  stationary  value.  Second,  these  cracks  propagate  in  a  transverse  direction 
along  the  weft  fibers  and  tend  to  cause  longitudinal  cracks  in  the  warp  direction.  The  third  and 
final  stage  is  that  of  longitudinal  crack  saturation  and  laminate  failure.  Apparently,  the  undulations 
are  the  “weak  links”  of  the  WR  laminates;  they  serve  as  local  stress  concentrators  due  to  fiber 
bending.  In  general,  WR  laminates  show  fatigue  resistance  superior  to  that  of  CSM  laminates  but 
inferior  to  the  non- woven  uni-directional  and  cross-ply  laminates  [Davis,  etal  (1964)].  Their 
primaiy  advantage  over  the  cross-ply  laminates  is  in  reducing  manufacturing  costs,  especially  when 
the  woven  fabric  is  used  in  combination  with  a  highly  efficient  fabrication  method  (e.g.,  the  resin 
transfer  molding  process). 

Obviously,  there  is  significant  difference  between  fatigue  failures  in  metals  and  in  FRP 
materials.  Fatigue  damage  in  metals  tends  to  localize  in  the  form  of  a  single  fatigue  macro-crack. 
As  it  extends,  a  crack  carries  a  severe  stress  concentration  at  its  tip  and  therefore  propagates  until 
it  reaches  a  region  of  low  stress  or  some  natural  obstacle  (e.g.,  a  bolt  hole).  Fatigue  damage  in 
FRP  laminates,  as  is  evident  from  the  previous  discussion,  is  much  more  diffuse.  Such  a  diffuse 
damage  property  inevitably  exhibits  local  reduction  of  elastic  moduli  which  attenuates  any  stress 
concentration  and  prohibits  damage  propagation.  (Experiments  with  specimens  containing  a  large 
rectangular  hole  with  saw-cut  comers,  [Smith  (1990)],  showed  that  following  one  million  load 
cycles,  the  only  evidence  of  fatigue  damage  were  whitened  patches  extending  approximately  1.25" 
from  the  hole  comers.)  In  general,  catastrophic  fatigue  failure  of  an  FRP  stmcture  by  through¬ 
thickness  crack  propagation  in  the  laminate  is  highly  unlikely  to  occur. 

Such  a  conclusion,  however,  does  not  suggest  that  stmctural  design  should  disregard 
fatigue  damage  accumulation  in  FRP  stmctures;  it  only  points  out  that  the  concern  for  fatigue 
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failure  should  concentrate  primarily  on  the  structural  joints.  Three  common  types  of  joints  used  in 
FRP  structures  are  bonded,  bolted,  and  bonded-bolted.  Fatigue  damage  in  any  of  these  joints  may 
manifest  in  different  ways  than  it  does  in  the  composite  laminate  itself  A  bonded  joint  may  fail  due 
to  a  crack  in  the  adhesive,  while  a  bolted  joint  may  exhibit  bearing  failure  of  the  bolt  hole.  In  any 
joint  configuration,  these  failures  may  be  supplemented  by  fatigue  damage  of  the  adherent  laminate 
plates.  Therefore,  fatigue  damage  processes  in  FRP  laminates  should  be  studied  in  the  context  of 
structural  joints.  This  conclusion  to  a  large  extent  motivates  the  selection  of  geometry  and 
configuration  of  the  test  specimens  described  later  in  the  experimental  section. 


MODELING  OF  FATIGUE  DAMAGE  IN  FRP  LAMINATES 

Modeling  of  fatigue  damage  has  two  major  aspects;  establishing  a  fatigue  failure  criterion 
and  a  damage  accumulation  rule.  In  the  traditional  approach  to  fatigue,  the  failure  criterion  relates 
fatigue  failure  (i.e.  fatigue  life)  to  the  number  of  applied  loading  cycles  to  cause  that  failure  (usually 
having  fixed  or  constant  amplitude).  An  analytical  expression  of  the  fatigue  failure  criterion  is  the 
S-N  curve,  which  plots  the  stress  level,  S,  of  the  constant-amplitude  cycles  against  the 
characteristic^  number  of  cycles  to  failure,  N,  at  this  stress  level.  In  modeling  fatigue  damage  with 
an  S-N  curve,  the  shape  of  the  curve  is  usually  fixed  by  an  empirical  parametric  expression;  these 
parameters  are  used  to  match  the  shape  to  a  particular  set  of  experimental  data.  Because  of  the 
large  scatter  of  experimental  fatigue  lives,  particularly  at  lower  stress  levels,  practical  measurement 


^  ‘Characteristic’  here  is  used  to  emphasize  that  fatigue  is  a  highly  stochastic  phenomenon, 
thereby  different  cycles  to  failure  may  be  associated  with  a  particular  stress  level.  In  a  most  general 
sense,  N  should  be  treated  as  random  number  and  assigned  probability  distribution.  Accordingly,  N 
is  usually  taken  as  the  ‘expected’  cycles  to  failure  but  can  be  the  cycles  to  failure  corresponding  to 
a  certain  probability  of  failure. 
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of  these  shape  parameters  requires  statistical  methods,  such  as  linear  and  nonlinear  regression 
analysis.  The  S-N  curve  alone  carmot  be  used  to  predict  fatigue  failures  under  variable-amplitude 
loadings.  A  damage  accumulation  rule  is  thus  postulated,  usually  by  prescribing  the  damage  done 
by  n  (n<N)  cycles  of  magnitude  S.  This  rule  is  often  formulated  independently  from  the  S-N  curve 
choice. 

Probably  a  more  sound  approach,  from  a  physical  point  of  view,  is  to  treat  damage 
accumulation  as  the  process  of  evolution  of  a  damage  parameter  under  the  loading  history.  This 
evolution  is  most  naturally  specified  in  terms  of  a  differential  equation  with  respect  to  the  loading 
cycles,  although  other  mathematical  expressions  may  also  prove  appropriate.  Examples  for  the 
damage  parameter  are  the  residual  stiffness  and  the  residual  strength  reviewed  later  in  the  report. 
In  association  with  such  an  approach,  the  failure  criterion  is  postulated  in  terms  of  a  critical 
(failure)  value  for  the  chosen  damage  parameter.  Integration  of  the  evolutionary  damage  equation 
can  be  done  for  both  constant  or  variable-amplitude  loadings.  For  the  constant  amplitude  case, 
analytical  integration  is  often  possible,  which  also  allows  a  relationship  between  the  critical  damage 
parameter  and  the  fatigue  life,  N,  to  be  established.  From  this  relationship,  an  S-N  curve  model 
can  also  be  derived,  which  puts  the  S-N  curve  approach  on  a  more  sound  basis. 

It  is  evident  from  the  previous  discussion,  that,  for  theoretical  as  well  as  practical  purposes, 
fatigue  damage  modeling  should  consider  both  constant  and  variable-amplitude  loading  situations. 
However,  we  prefer  to  discuss  S-N  curves,  damage  parameters,  evolutionary  and  failure  equations, 
and  all  other  issues  with  constant-amplitude  damage  accumulation  in  FRP  laminates  in  this  section 
of  the  report.  The  next  section  of  this  report  is  devoted  to  the  specific  issues  associated  with 
fatigue  damage  modeling  under  variable-amplitude  loadings. 

Among  all  structural  materials,  fatigue  of  metals  has  been  studied  most  extensively,  and  it 
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has  been  found  that,  almost  exclusively,  a  metal’s  S-N  curve  is  best  expressed  in  the  form  of  a 
power  law  as 


where  m  and  K  are  known  as  the  fatigue  exponent  and  the  fatigue  coefficient.  These  are  material 
parameters,  determined  usually  from  a  linear  regression  fit  of  log  (S)  versus  log  (N)  fi'om 
experimental  data.  Another  power-law  relationship  established  in  metal  fatigue  is  the  Manson- 
Coffin  relationship,  usually  expressed  as 


Bap  BfO^ )  (2) 

in  which  Sap  is  the  applied  plastic  strain  amplitude,  and  e/  and  c  are  fatigue  coefficients.  The 

Manson-Coffin  equation  is  usually  employed  in  the  context  of  local  approaches  to  fatigue  damage 
accumulation;  only  in  such  a  situation  can  plastic  strain  magnitude  be  quantified  from  local  material 
response  due  to  stress  concentration  and/or  residual  stresses. 

Because  it  is  difficult  to  define  plastic  strain  for  composite  materials,  the  Manson-CoflHn 
relation  has  not  found  any  application.  However,  due  to  its  universality,  the  power  law  S-N  curve 
is  often  used  to  characterize  fatigue  failures,  especially  when  little  information  on  the  fatigue 
behavior  of  a  particular  material  exists.  In  composite  fatigue,  this  relationship  is  usually  associated 
with  the  names  of  Hahn  and  Kim  [1976]  who  preferred  to  express  it  as 

sN^=^c  (3) 
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In  this  equation,  d  and  c  are  model  (material)  parameters,  and  s  is  the  ratio  of  the  applied 
stress,  S,  and  the  material’s  ultimate  strength.  Suit,  i.e., 

S 

s  =  — 

S  uh  ^4) 

This  ratio  has  found  a  considerable  application  as  a  measure  of  the  loading  severity  in 
fatigue  studies.  Hahn  and  Kim  also  consider  the  semi-logarithmic  S-N  curve  model  in  which  the 
applied  stress  level  is  related  to  log(N).  This  is  expressed  as 

s  =  k\og(N)  +  b 

where  k  and  b  are  material  parameters  and  s  is  given  by  Equation  (4). 

These  two  models  are  compared  in  [Hahn  and  Kim  (1976)]  for  experimental  data  obtained 
from  [0/±  45/90]  glass/epoxy  laminates  for  up  to  10^  cycles.  It  was  concluded  that,  overall. 
Equation  (3)  gives  a  better  fit  than  Equation  (5)  when  the  two  models  are  compared  over  the  entire 
range  of  fatigue  lives  obtained,  although  in  the  low-cycle  fatigue  range  Equation  (5)  can  provide 
superior  results. 

Another  simple  analytical  S-N  curve  model  is  that  proposed  by  Hwang  and  Han  [1986a, 
1986b].  The  main  difference  between  this  model  and  those  previously  specified  is  that  it  represents 
in  inverse  power  law,  i.e., 

N  =  [B(l-s)f^ 
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in  which  B  and  C  are  model  parameters.  Equation  (6)  is  derived  as  a  particular  case  of  the  fatigue 
modulus  concept  combined  with  a  strain  failure  criterion;  both  are  discussed  below  in  more  detail. 

Based  on  the  fatigue  data  reported  by  Hwang  and  Han  [1986b]  for  aerospace  types  of  FRP 
laminates,  and  on  that  reported  by  Howson,  etal.  [1992],  and  Owen  [1968]  for  marine  types  of 
composites.  Read  and  Shenoi  [1995]  performed  a  comparison  of  the  models  specified  by  Equations 
(3),  (5),  and  (6).  They  plotted  all  test  data  on  a  semi-logarithmic  scale.  (On  such  a  scale  Equation 
(5)  plots  as  a  straight  line.  Equation  (3)  plots  as  a  bottom  convex,  and  Equation  (6)  plots  as  a 
bottom  concave.)  The  data  from  Hwang  and  Han  [1986b],  obtained  from  glass/epoxy  laminates,  is 
of  typical  low-cycle  fatigue  type  (less  than  10^  cycles).  For  these  data,  the  best  fit  was  provided  by 
Equation  (6),  but  it  was  only  marginally  better  than  Equation  (5).  Howson,  etal,  [1992]  reports 
results  from  two  types  of  composites:  one  made  of  glass  fibers  in  polyester,  and  another  made  of 
hybrid  glass-aramid  fibers.  The  longest  fatigue  lives  obtained  are  above  10^  but  less  than  10^ 
cycles.  For  both  composites  it  was  concluded  that  Equation  (5)  provides  the  best  overall  fit. 
However,  it  could  be  easily  observed  that  the  fit  was  particularly  poor  for  low-  and  high-cycle 
fatigue.  For  the  low-cycle  fatigue  data  alone  (less  than  20  cycles)  Equation  (3)  gave  the  best  fit. 
Owen  and  Howe  [1968]  report  high-cycle  fatigue  results  from  glass/polyester  CSM  composites. 
The  best  fit  was  again  obtained  with  Equation  (6). 

Although  fatigue  data  for  typical  marine  FRP  laminates  are  relatively  limited  compared  to 
aerospace  laminates,  it  is  evident  fi-om  the  quoted  investigations  that  the  fatigue  behavior  of  the 
variety  of  FRP  types  encountered  in  marine  applications  may  not  be  characterized  by  a  single 
anal5^ical  form  of  the  S-N  curve.  Moreover,  which  S-N  curve  model  is  to  be  used  may  have  to  be 
decided  for  each  particular  material.  Therefore,  it  is  highly  desirable  to  have  a  more  general  fatigue 
damage  model,  based  on  the  micro  mechanisms  of  composite  damage,  from  which  particular 
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damage  models  can  be  formulated  to  fit  particular  applications.  In  FRP  fatigue,  theoretical 
background  for  such  a  model  is  provided  by  three  somewhat  related  concepts:  stiffness 
degradation,  fatigue  modulus  degradation,  and  strength  degradation. 

The  concept  of  stiffness  degradation  stems  directly  from  the  diffuse  character  of  the  fatigue 
damage  in  composite  laminates,  as  discussed  previously.  On  the  macroscopic  scale,  the  effects  of 
the  numerous  matrix  cracks,  fiber-matrix  debonds,  fiber  breaks,  and  ply  delaminations  result  in 
reduction  of  stififiiess.  The  reduction  of  stiffness  can  therefore  be  used  as  a  measure  of  the  damage 
developed.  The  residual  stiffness  approach  was  proposed  by  Yang,  et  al.  [1989,  1990].  The  basic 
assumption  of  the  approach  is  that  the  rate  of  degradation  of  residual  stiffhess  is  a  power  function 
of  the  number  of  cycles.  This  can  be  expressed  as 


dEn 

dn 


=  -EoQvn' 


(7) 


in  which  Eo  and  £„  are  the  initial  stiffhess  and  the  stififiiess  after  n  applied  cycles,  respectively.  Q 
and  V  are  model  parameters  which  account  for  the  influence  of  applied  stress  level,  loading 
frequency,  and  environment.  Integration  of  Equation  (7)  renders: 


E„  =  EoO-Qn) 

Based  on  the  experimental  data,  it  was  suggested  that  Q  depends  on  v  linearly,  i.e. 


(8) 


Q-ai^  a2'V 


(9) 
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and  that,  v ,  on  its  part  depends  linearly  on  the  applied  stress  level,  S,  i.e. 

\  =  a3  +  BS  (20) 

By  using  Equations  (9)  and  (10),  Equation  (8)  can  be  expressed  as 


En""  Eo[l-(  ai  +  a2a3  +  ] 


(11) 


Thus  four  parameters,  ai,  a2  ,a3,  and  B,  independent  of  the  applied  stress  level  S,  completely  define 
the  above  residual  stiffness  model. 

If  a  failure  criterion  is  postulated,  an  S-N  curve  expression  based  on  Equation  (11)  can  be 
developed.  In  a  later  study,  Lee,  et  al.,  [1996]  postulated  that  failure  takes  place  when  critical 
strain  is  reached.  This  study  actually  deals  with  variable-amplitude  damage  accumulation  but,  as 
explained  in  the  next  section,  the  specifics  of  the  approach  allows  it  to  be  applied  to  constant- 
amplitude  situations  as  well.  Accordingly,  the  strain  at  the  failure  cycle  N  is  defined  somewhat 
empirically  as 


^  Sun 
ew-  — +  ao— 
Eo  Eo 


(12) 


From  Equation  (11),  the  fatigue  life,  N,  can  be  expressed  as  a  function  of  the  stiffhess  at  failure, 
En,  as 

(13) 


N  = 


1  ~  En  ^  Eo 

ai  +  02  03'^  a2B  S 
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Considering  that  ^ En  and  taking  into  account  Equation  (12),  the  fatigue  life  can  be 


expressed  as 


N  = 


~\l/(  S) 


1-- 


‘S'+flo  Suit 


ai  ^  at  02'^  a2 


BS 


(14) 


which  is  the  sought  after  S-N  curve  formula.  Evidently,  it  provides  for  five  independent  model 
parameters,  ao,  ai,  a2  ,a3,  and  B,  which  are  to  be  determined  from  experimental  data.  Obviously, 
this  is  a  much  more  involved  task  than  determining  the  two  coefficients  of  the  power  law  S-N 
curve  (3).  To  this  end,  according  to  the  procedure  suggested  by  Yang,  etal,  [1990],  residual 
stiffhess  is  measured  during  a  constant  amplitude  test  at  a  particular  stress  level:  Ei,  E2,  ...,  Ek.  at 
cycles  ni,  n2,  ...,  Uk.  respectively.  By  linear  regression  of  the  residual  stiffnesses,  Ek,  with  respect  to 
the  loading  cycles,  nk,  the  model  parameters  Q  and  v  can  be  determined  as  evident  from  Equation 
(8).  By  subsequent  regressions  with  respect  to  different  loading  levels,  S,  the  parameters  ai,  a2  ,a3, 
and  B  can  also  be  found  as  suggested  by  Equations(9)  and  (10). 

It  is  worth  noting  that  in  the  approach  described  by  Yang,  et  al.,  [1989,  1990],  the  residual 
stiffhess  approach  is  applied  in  a  probabilistic  form.  In  particular,  the  initial  stiffhess,  Eo,  and  the 
two  model  parameters,  Q,  v,  and  B  are  considered  to  be  random  variables.  At  the  same  time,  ao, 
ai,  a2,  and  as  are  still  treated  as  deterministic.  The  experiments  indicated  that  the  distributions  of 
both  Eo  and  B  are  reasonably  well  modeled  as  log-normal.  Thus,  the  probability  distribution  of  Eo 
is  fully  characterized  by  the  mean  and  the  standard  deviation  of  ln( Eo),  and  Oo;  and  that  of  B, 
by  the  mean  and  the  standard  deviation  of  \n(B),  |4b  and  ob.  Accordingly,  the  probability 
distribution  of  the  residual  stiffhess  at  the  nth  cycle,  En,  depends  on  three  deterministic  parameters, 
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ai,  a2,  and  as,  and  on  the  probability  distributions  of  two  random  variables,  Eo  and  B.  This  brings 
to  seven  the  total  number  of  independent  parameters  of  the  model  (ai,  a2,  as,  po,  cto,  Pb,  and  CTb)  if 
the  interest  is  on  the  probability  distribution  of  the  residual  stiffiiess  only. 

With  this  objective  Yang,  et  al  [1990],  applied  the  residual  stififtiess  to  analyze  experimental 
data  from  constant-amplitude  fatigue  tests  of  graphite/epoxy  [0°,  ±45°,  90°]  laminates.  They 
compared  the  experimentally  obtained  distribution  of  the  residual  stiffness  with  the  theoretical 
obtained  through  Equation  (11).  Fairly  good  correlation  between  the  experimental  results  and  the 
theoretical  predictions  was  reported. 

If  the  objective  is  to  compute  the  parameters  of  the  S-N  curve  (14),  the  initial  stiffness  Eo 
does  not  have  to  be  estimated,  but  the  static  ultimate  strength.  Suit,  must  be  considered.  In  [Halpin, 
etal  (1973)],  Suit  was  modeled  as  a  random  variable  whose  probability  distribution  is  a  two- 
parameter  Weibull,  i.e. 


PsJS)  = 


a 

a-l 

fs] 

a  ’ 

Uj 

exp 

(15) 


If  the  theorem  of  total  probability  is  applied  to  Equation  (14)  along  with  the  expressions  for 
the  probability  distributions  of  B  and  Suit,,  the  probability  distribution  of  N  and  any  other  statistical 
parameter  (e.g.,  its  expected  value)  can  be  found.  The  total  number  of  material  parameters  that 
need  to  be  estimated  from  testing  is  eight  (ao,  ai,  a2,  as,  Pb,  CTb,  ct  and  P).  Note  that  only  the  first 
six  of  them  are  computed  from  the  fatigue  tests,  while  the  distribution  parameters  of  the  ultimate 
static  strength,  ot  and  p,  are  obtained  from  the  static  tests.  In  general,  Yang  s  residual  stiffness 
approach  is  quite  complex  and  experimentally  expensive,  and  therefore  should  be  considered  only  if 
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other  more  simple  approaches  fail  to  provide  satisfactory  modeling  of  fatigue  damage. 

Based  on  the  notion  of  degrading  stiffness,  Hwang  and  Han  [1986a  and  1986b]  proposed 
the  concept  of  fatigue  modulus  that  is  different  from  the  stiffness  (elastic  modulus).  Accordingly, 
fatigue  modulus  at  the  nth  cycle,  Fn,  is  defined  as 


F„  = 


S 


Sn 


(16) 


where  e„  is  the  resultant  strain  at  the  maximum  load  of  the  cycle.  The  basic  assumption  is  again 
that  the  modulus  degradation  rate  is  a  function  of  the  number  of  fatigue  cycles,  which  can  be 
expressed  in  a  form  similar  to  Equation  (7)  with  stiffness  replaced  by  fatigue  modulus  [Hwang, 
etal,  (1986b)].  More  recently,  it  was  suggested  that  this  degradation  rate  is  a  function  of  both  the 
number  of  fatigue  cycles  and  the  current  value  of  the  fatigue  modulus  [Hwang,  etal,  (1995)].  By 
following  this  more  general  assumption,  the  fatigue  modulus  degradation  rate  can  be  expressed  as 


in  which  A,  B,  and  C  are  model 
this  general  equation  renders 


dn  B  F^'‘ 


(17) 


parameters  to  be  obtained  from  experimental  data.  Integration  of 


Fl-Fl  =  -Ari 


Thus,  the  number  of  cycles  to  failure,  N,  can  be  expressed  as 


(18) 
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(19) 


N  =  [(F^o-F})/Af^ 

in  which  Ff  is  the  fatigue  modulus  at  failure.  Note  that  the  above  equation  simply  associates  the 
number  of  cycles  with  the  fatigue  modulus  at  failure,  but  does  not  relate  either  of  them  to  the 
magnitude  of  the  stress  cycle.  A  way  out  of  this  difficulty  is  offered  again  by  the  failure  criterion. 
Similarly  to  Equation  (12),  it  is  possible  to  associate  failure  with  the  advent  of  failure  strain.  In  a 
more  simple  form,  the  failure  criterion  can  postulate  that  fatigue  failure  occurs  at  the  same  strain  as 
the  static  failure  of  undamaged  material,  i.e., 

Po  (20) 

Another  approach,  chosen  by  Hwang,  et  al.,  [1986b],  is  to  model  failure  indirectly,  by  addition  of 
another  model  parameter.  To  this  end,  they  introduce  a  reference  fatigue  modulus,  Fr,  such  that 

El=J-=s  p±=p 

Fr  Suit  Fr  (21) 

It  can  be  shown  that  the  p-parameter  in  the  above  equation  has  a  simple  physical  meaning:  it  is  the 
ratio  of  the  strain  at  fatigue  failure,  £/,  to  the  strain  at  static  failure  of  undamaged  material,  . 

Accordingly,  the  case  p=l  corresponds  to  the  strain  failure  criterion,  Equation  (20).  The 
experimental  data  presented  by  Hwang,  etal,  [1995]  demonstrates  that,  in  all  cases,  p  is  quite 
close  to  unity. 

The  fatigue  life.  Equation  (19),  can  now  be  expressed  as 


(22) 


N-[M(p^-sVf'' 

where  M  =  f^/ A.  Although  developed  from  the  fatigue  modulus  concept,  the  above  equation 
has  no  explicit  use  of  fatigue  modulus  parameters.  Instead,  it  contains  four  model  parameters  (p, 
M,  B,  C)  that  should  be  determined  from  the  constant-amplitude  experimental  data. 

Again,  Equation  (22)  provides  a  generalized  S-N  curve  model  for  fatigue  life  under 
constant-amplitude  loading.  Hwang,  etal,  [1995]  refer  to  it  as  the  modified  fatigue  life  prediction 
equation  (MFLPE).  Note  that  for  p=l  (strain  failure  criterion)  and  B=1  it  results  in  the  inverse 
power  law  S-N  curve,  Equation  (6).  Other  S-N  curve  models  can  also  be  derived  from  the  fatigue 
modulus  concept,  although  not  directly  from  Equation  (17)  or  Equation  (22).  Consider  the 
following  particular  forms  of  the  general  fatigue  modulus  degradation; 

^  =  -A-  (a)  ^  =  -—  (b) 

dn  n  dn  n  (23) 

Integrating  the  above  equations  (this  time  taking  n  to  vary  from  1  to  N)  and  taking  into  account 
Equation  (21),  renders  for  the  fatigue  cycles  at  failure 


N  = 


(zl)  N  =  exp  [m(p®  -  s®)]  (b) 


(24) 


respectively.  The  first  of  these  equations  is  a  form  of  the  power  law  S-N  curve.  Equation  (3), 
while  the  second  is  a  form  of  the  semi-logarithmic  S-N  curve.  Equation  (5).  Thus,  Hwang’s 


19 


formula  indeed  provides  quite  a  general  framework  for  development  of  different  fatigue  damage 
models. 

Fatigue  damage  within  a  composite  laminate  can  not  only  degrade  the  elastic  moduli,  but 
can  also  lower  its  static  strength.  Based  on  this  observation,  the  residual  strength  was  proposed  by 
Halpin,  etal,  [1973]  as  the  physical  parameter  describing  damage  in  composites.  Much  work  on 
the  residual  strength  approach  has  been  done  by  Yang,  etal,  [1980,  1981  and  1983].  In  particular 
they  suggested  that  an  equation  for  the  rate  of  change  of  the  residual  strength  Sr  during  the  nth 
loading  cycle  can  be  expressed  as 


dSr_ 

dn  m  (25) 

in  which  A  and  m  are  model  parameters.  It  was  also  suggested  that  A  depends  on  the  applied 
stress  at  the  nth  cycle,  S,  and  the  initial  static  strength.  Suit.  By  integrating  this  equation, 
expressions  similar  to  those  developed  for  the  residual  stiffness  and  the  fatigue  modulus  concept 
can  be  derived. 

Literature  reveals  that  the  residual  strength  concept  was  proposed  earlier  than  those  of 
residual  stiffness  and  fatigue  modulus.  It  seems,  though,  that  the  residual  strength  concept  has  lost 
popularity  among  fatigue  analysts  at  the  expense  of  the  residual  stiffiiess  approach  (and  the  related 
fatigue  modulus  approach).  The  reasons  are  pointed  out  by  Hashin  [1985].  Clearly,  because  the 
residual  stiffness  is  a  non-destructive  parameter,  it  can  not  only  be  calculated  from  specimen 
response  without  even  interrupting  the  test  but  can  also  be  measured  on  a  structure  of  interest. 
Moreover,  residual  stiffhess  can  in  principle  be  determined  analytically  from  a  given  crack 
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distribution.  It  has  also  been  demonstrated  that  residual  stifi&iess  exhibits  much  less  statistical 


scatter  than  residual  strength.  Finally,  it  may  be  noted  that  within  Kachanov’s  concept  of 
‘continuum  damage’,  the  residual  strength  approach  would  be  equivalent  to  the  residual  stiffiiess 
approach  as  long  as  the  material  does  not  exhibit  significant  ductility  (which  is  usually  the  case  with 
FRP  composites).  Thus,  it  is  highly  unlikely  if  any  practical  model  for  fatigue  damage 
accumulation  in  FRP  laminates  can  be  developed  exclusively  on  the  residual  strength  concept. 


MODELING  DAMAGE  ACCUMULATION  IN  FRP  LAMINATES 

Cumulative  damage  theory  is  probably  the  only  approach  that  can  predict  the  fatigue  life  of 
a  composite  laminate  under  a  multi-level  stress  history  because  application  of  fracture  mechanics 
techniques  is  highly  questionable  due  to  the  diffuse  nature  of  fatigue  damage  in  composites.  In 
general,  the  cumulative  damage  approach  postulates  that  a  damage  function,  D,  exist  such  that 


D  = 


0  material  never  subjected  to  loading 
1  at  failure 


(26) 


In  addition  to  the  above  requirements,  damage,  D,  is  postulated  to  be  a  monotonically  increasing 
function  of  the  loading. 

For  modeling  damage  accumulation  due  to  fatigue,  the  main  objective  is  to  analytically 
express  damage,  D,  caused  by  a  segment  of  n  constant-amplitude  cycles  with  magnitude  S  given 
that,  the  applied  cycles,  n,  is  less  than  the  constant-amplitude  fatigue  life,  N.  Following  Hwang, 
etal,  [1986a],  for  a  composite  material  it  can  be  most  generally  expressed  as 
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D^D(Do.  n,  S,  f  T,  M,...) 


(27) 


where  Do  is  initial  damage  existing  prior  to  the  application  of  the  considered  loading  segment,  n  is 
the  number  of  fatigue  cycles,  S  is  the  applied  stress  level,  f  is  the  loading  frequency,  T  is  the 
temperature,  and  M  is  the  moisture  content.  The  effect  of  the  last  three  parameters  on  fatigue 
damage  will  not  be  considered  in  this  section. 

Considering  the  combined  influence  of  the  other  three  parameters.  Do,  n,  and  S,  on  the 
fatigue  damage  accumulation,  damage  is  assumed  to  be  caused  by  a  multi-level  stress  history. 
Therefore,  not  only  does  the  damage  contribution  from  each  loading  cycle  need  to  be  added,  but 
the  influence  of  the  loading  of  each  cycle  on  the  damage  caused  by  subsequent  cycles  needs  to  be 
identified,  as  well.  This  is  quite  a  complex  situation,  which  is  rarely  employed  in  practice.  The  first 
simplification  is  done  by  neglecting  the  effect  of  Do;  that  is,  by  assuming  that  damage  caused  by  the 
sequence  of  n  cycles  with  magnitude  S  is  always  the  same  regardless  of  where  this  sequence  occurs 
in  the  loading  history.  Such  damage  accumulation  rules  have  been  considered  in  the  context  of 
fatigue  of  FRP  laminates  and  are  discussed  below. 

The  simplest  damage  accumulation  law,  which  is  commonly  known  as  the  Palmgren-Miner 
rule,  [Palmgren  (1924),  Miner  (1945)],  assumes  that  damage  contribution  from  each  cycle  of  the 
loading  history  is  independent  from  the  other  cycles.  Therefore,  the  damage  inflicted  by  n  stress 
cycles  with  magnitude  of  S  is  simply 

D  =  — 

N  (28) 

where  N  denotes  the  cycles  to  failure  at  S  from  the  constant-amplitude  S-N  curve.  For  all  stress 
levels  this  damage  rule  yields 
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(29) 


Hi 

Ni 


where  iii  is  the  number  of  cycles  having  amplitude  Si. 

In  many  situations,  the  rate  of  fatigue  damage  accumulation  in  composites  may  deviate 
significantly  from  that  assumed  in  Equation  (28).  For  example,  Owen  and  Howe  [1968]  report 
that 
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(30) 


fits  much  better  constant-amplitude  data  they  obtained  from  CSM  composites.  Such  an  equation, 
however,  has  no  particular  benefit  over  the  linear  damage  accumulation  rule,  Equation  (28). 
Making  damage  a  function  of  the  ratio  of  n/N  only  implies  that  it  is  always  the  same  regardless  of 
the  cycle’s  position  in  the  loading  sequence.  Therefore,  there  is  always  a  measure  of  damage 
which  accumulates  linearly  under  constant-amplitude  loadings. 

In  connection  with  Equation  (1),  the  Palmgren-Miner  rule.  Equation  (29),  is  used  almost 
exclusively  in  metal  fatigue.  By  treating  a  multi-level  stress  time-history  as  a  random  process. 
Equations  (1)  and  (29)  allow  a  simple  stochastic  approach  to  be  used  to  estimate  cumulative 
damage.  This  approach  is  applicable  when  the  loading  can  be  described  as  a  zero-mean,  Gaussian, 
and  narrowband  stochastic  process.  In  such  a  case,  the  amplitudes  of  the  fatigue  cycles  possess  a 
Rayleigh  probability  distribution,  which  can  be  expressed  as 


Ps(S)  = 


O  rms 


exp 


2S 


,2 

rms  j 


(31) 
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In  the  above  formula  Snm  denotes  the  root-mean-square  value  of  the  stochastic  stress  history.  By 
using  the  Palmgren-Miner  hypothesis,  Equation  (28),  and  the  power-law  S-N  curve.  Equation  (1), 
the  expected  fatigue  damage  E[D]  and  the  expected  cycles  to  failure  E[N]  are  expressed  as 


E[D]  = 
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E[N] 


K  ^  ^  k{s\ 
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The  integration  of  the  above  equation  renders 
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where  r(*)  is  the  gamma  function.  Equation  (33)  is  known  as  the  “Rayleigh  approximation” 
equation  [Miles  (1954)]. 

An  expression  similar  to  Equation  (33)  can  be  derived  if  the  Palmgren-Miner  rule  is  applied 
in  conjunction  with  the  semi-logarithmic  S-N  curve  model.  Equation  (5).  To  this  end,  assume  first 
that  the  semi-logarithmic  model  is  expressed  as 


\og(N)  =  b-aS 


(34) 


in  which  b  and  a  are  the  parameters  of  the  model  obtained  from  a  linear  regression  of  the 
experimental  data.  This  equation  can  be  further  rearranged  as 
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where 


(35) 


N  =  jQb-aS  ^  ^b-aS 


b  =  b\n(10)  a  =  a\n(10) 


(36) 


The  expected  damage  can  then  be  expressed  as 
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Performing  the  integration  in  Equation  (37)  gives 


E[D]  = 
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(38) 


in  which  0(*)  denotes  the  Gaussian  cumulative  distribution  function. 

Damage  accumulation  rules  in  which  D  cannot  be  expressed  as  a  function  of  n/N  have  also 
been  applied  to  composite  materials.  A  compilation  of  such  models  can  be  found  in  Hwang,  etal, 
[1986a].  The  damage  function  in  such  a  case  can  be  expressed  as  D(n,S),  as  suggested  in  Equation 
(29),  or  alternatively  as  D(n,N)  because  S  is  always  a  single- valued  ftmction  of  N.  Derivation  of 
such  models  can  be  based  on  pure  empirical  considerations  from  particular  experimental  data  or  on 
the  grounds  of  a  more  physical  damage  hypothesis  such  as  the  residual  strength,  residual  stififhess, 
or  fatigue  modulus  concepts  reviewed  in  the  previous  section. 

An  example  of  the  empirical  approach  is  the  damage  accumulation  rule  suggested  by 
Hashin  [1985].  Based  on  observations  that  fatigue  damage  varies  nearly  logarithmically  with  the 
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number  of  cycles,  he  suggests  that  damage  accumulated  after  n  constant-amplitude  cycles  with 
magnitude  S  is 


\og(n/Ne) 

\og(N/Nj 


(39) 


in  which  Ne  is  a  model  parameter. 

An  example  of  the  physical  approach  is  a  damage  accumulation  model  based  on  the  residual 
strength  concept  also  proposed  by  Hashin  [1985],  The  model  is  expressed  as 


1-S  (40) 

in  which  Sr  is  the  ratio  of  the  residual  strength,  Sr,  to  the  ultimate  static  strength,  Suu.  The 
relationship  between  Sr  and  the  applied  number  of  cycles,  n,  is  provided  by  Equation  (25). 

Another  damage  model  of  this  type  is  proposed  by  Hwang  and  Han  [1986b].  This  model  is 
based  on  Equations  (23b)  and  (24b)  and  is  expressed  as 


D  = 


f  ^  1 

(41) 


The  model  parameters  B  and  C  in  the  above  equation  are  found  from  the  inverse  power  law  S-N 
curve.  Equation  (6). 

Many  possibilities  for  development  of  nonlinear  damage  accumulation  laws  exist  when  the 
process  is  based  on  the  degradation  models  discussed  in  the  previous  section.  Hwang,  et  al.. 
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[1986b,  1995]  demonstrate  this  in  the  context  of  the  fatigue  modulus  concept.  In  Hwang,  etal, 
[1995],  they  considered  four  different  types  of  damage  accumulation.  According  to  these  laws,  the 
accumulated  fatigue  damage  after  n  cycles  with  amplitude  S  is  taken  as 

B  B 

£,  =  EozIjl(1)  d  =  ^(II)  D=^^^(III)  D  =  ^^^(IV) 

Fo-Ff  Ef  ecEo  Fo-Ff  ^42) 

In  the  above  equation,  (Fo,  ),  (Fn,  e„ ),  and  (Ff,  e/ )  are  the  fatigue  modulus  and  the  maximum 

strain  at  the  beginning,  after  n  cycles,  and  at  failure  under  constant-amplitude  loading,  respectively. 
Application  of  the  above  damage  accumulation  hypotheses  to  different  degradation  models  derived 
as  particular  cases  of  Equations  (19)  and  (24)  allowed  these  damage  models  to  be  expressed  in 
terms  of  the  p,  M,  A,  B,  and  C  parameters  defined  in  Equations  (19)  and  (24).  In  particular. 
Equation  (19)  expresses  the  damage  models  of  Equations  (42)  as 
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Use  of  nonlinear  damage  models  such  as  Equations  (39)-(43)  in  multi-level  stress  histories 
cannot  be  done  by  simply  adding  up  damage  due  to  cycles  with  different  stress  levels.  To  perform 
the  summation,  a  procedure  based  on  the  concept  of  ‘equivalent  damage’  is  usually  applied.  To 
explain  the  procedure,  assume  that  damage  is  specified  in  the  form  D(n,S)  and  that  a  loading 
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history  consisting  of  ni  cycles  of  magnitude  Si,  n2  cycles  of  magnitude  S2,  ns  cycles  of  magnitude 
S3,  and  so  on,  is  applied.  After  the  application  of  ni  cycles  of  magnitude  Si  the  accumulated 
damage  is 

Di=  D(ni,Si)  Di<l  (44) 

The  damage-equivalent  cycles  n2e  of  magnitude  S2  which  cause  an  amount  of  damage  equal  to  the 
accumulated  damage,  Di,  is  found  by  solving  the  following  equation: 


Dl-  D(n2e’S2)  (45) 

The  accumulated  damage  after  the  application  of  the  second  sequence  of  nj  cycles  of  magnitude  S2 
is  now  found  as 

D2  =  D( n2  +  n2e  .  S2)  D2<  1  (46) 

The  procedure  is  repeated  until  failure  is  achieved  (Dk=l)  or  the  loading  sequence  is  exhausted. 

Such  a  procedure,  although  conceptually  simple,  may  actually  be  quite  difficult  to  apply  in 
practice.  First,  during  each  step  of  obtaining  damage-equivalent  cycles.  Equation  (45),  a  nonlinear 
equation  may  have  to  be  solved.  This  can  make  the  procedure  quite  expensive  computationally, 
given  that  maritime  loadings  of  interest  may  contain  millions  of  cycles  before  final  failure  is 
achieved.  A  much  more  serious  problem,  however,  is  represented  by  the  broadband  loadings.  A 
load  cycle  identification  scheme  such  as  the  rainflow  algorithm  may  be  used  to  determine  the 
individual  fatigue  cycles  of  a  broadband  loading  history.  The  lack  of  order  in  the  fatigue  cycles 
obtained  by  such  a  scheme,  however,  makes  the  rainflow  algorithm  practically  inapplicable  in 
connection  with  damage  models  like  those  considered  above.  In  other  words,  the  damage 
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accumulation  theory  does  not  possess  analytical  means  to  compute  accumulated  fatigue  damage 
from  nonlinear  models  and  broadband  loading  histories. 

A  possible  way  out  of  this  difficulty  may  be  found  by  systematic  application  of  probabilistic 
approaches  (random  vibration  methods  in  particular)  to  evolutionary  damage  equations  such  as  the 
degradation  models  specified  by  Equations  (7),  (17)  and  (25).  However,  application  of 
probabilistic  approaches  to  the  problem  of  composite  fatigue  under  variable-amplitude  loadings  are 
quite  rare.  One  of  the  few  studies  is  that  of  Lee,  et  al.,  [1996]  who  applied  Yang’s  residual 
stiffness  model.  Equations  (7)-(15),  to  variable-amplitude  fatigue  damage  accumulation  under 
block-type  spectrum  loadings.  Specifying  the  fatigue  loading  in  terms  of  a  block-type  spectrum  is 
quite  common  in  aerospace  applications,  however,  is  rarely  used  in  marine  applications.  Rather 
than  evaluating  the  parameters  of  the  model  from  constant-amplitude  loadings,  they  applied 
Equations  (11)  and  (14)  directly  to  fatigue  data  obtained  under  the  block  loadings  without  any 
modification.  In  this  approach,  however,  n  refers  not  to  the  nth  loading  cycle  but  to  the  nth 
loading  block  in  the  spectrum.  Accordingly,  N  designates  the  block  loading  number  at  which 
failure  occurs.  The  stress  variable  S  in  Equation  (11)  is  postulated  to  be  equal  to  the  maximum 
stress  in  the  block.  The  other  assumptions  of  the  model,  particularly  those  on  the  probability 
distributions  of  Eo,  Suit,  and  B  are  preserved.  To  validate  this  approach,  experiments  with 
[0°/907  ±  45°]  graphite/epoxy  laminates  were  carried  out.  Four  block  loading  sets  with  different 
stress  levels  were  generated.  In  predicting  the  statistical  distribution  of  the  fatigue  life  for  any  one 
of  the  data  sets,  the  remaining  three  data  sets  were  considered  as  a  basic  data  from  which  the 
model  parameters  were  estimated.  The  theoretical  predictions  were  then  correlated  with  the 
experimental  results  of  that  data  set.  Fairly  good  agreement  between  the  predictions  and 
experiments  was  found.  These  positive  results,  however,  should  be  taken  as  evidence  for  the 
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applicability  of  the  residual  stiffness  model  to  composite  fatigue.  At  the  same  time,  it  seems 
dubious  that  the  method  described  above  can  be  applied  directly  to  applications  which  do  not 
specify  the  fatigue  loading  as  block-type  spectra. 

Unfortunately,  experimental  data  on  cumulative  damage  is  mostly  related  to  composites  in 
aerospace  applications.  In  fact,  the  majority  of  the  data  is  for  two-stage  loadings  and  low-cycle 
fatigue.  For  such  data  Read  and  Shenoi  [1995]  compare  some  of  the  above  damage  accumulation 
models.  In  particular,  they  considered  the  Palmgren-Miner  rule,  Kashin’s  damage  model. 
Equation  (39),  (which  for  two  stage-loadings  and  Ne==l  is  the  same  as  a  method  presented  by 
Kashin  and  Rotem  [1978]),  Kashin’s  residual  strength  model.  Equation  (40),  and  the  fatigue 
modulus  model  of  Kwang  and  Kan  ,  Equation  (41).  They  conclude  that,  in  most  cases,  the 
predictions  made  using  Kashin  and  Rotem’ s  method,  are  closest  to  experimental  data,  while  the 
approach  of  Kwang  and  Kan  tends  to  give  the  most  inaccurate  predictions.  In  many  instances,  the 
predictions  of  the  Palmgren-Miner  approach  were  close  to  those  of  Kashin  and  Rotem’ s  method. 
Overall,  the  linear  damage  accumulation  rule  made  predictions  which  were  on  the  non-conservative 
side. 

Based  on  the  studies  reviewed,  the  only  sure  conclusion  that  can  be  made  is  that  sufficient 
experimental  data  pertinent  to  fatigue  of  FRP  laminates  under  multi-level  stress  histories  simply  do 
not  exist.  Therefore,  it  is  too  early  at  this  stage  to  draw  conclusions  about  which  damage 
accumulation  approaches  will  render  the  most  accurate  results  in  marine/naval  applications— all 
options  should  be  considered  in  the  study. 
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EXPERIMENTAL  EVALUATION  OF  COMPOSITE  LAMINATE  AND  JOINT  DETAILS 


The  composite  test  program  discussed  in  this  report  was  set  up  to  provide  a  realistic 
insight  into  the  static  and  fatigue  behavior  and  damage  accumulation  of  FRP  laminates  and  joints. 
To  study  the  fatigue  damage  accumulation,  several  different  specimen  configurations  were 
designed  and  fabricated  for  characterization  under  constant  and  variable  amplitude  fatigue 
loadings.  The  objective  of  the  tests  on  the  base  specimens  was  to  investigate  the  fatigue  strength  of 
a  bi-directional  base  laminate  under  four  different  axial  load  scenarios.  Test  scenarios  included  two 
loading  orientations,  0/90°  and  ±  45°,  with  respect  to  the  laminate  lay-up,  and  two  different  cases 
of  load  transfer.  One  case  of  load  transfer,  referred  to  as  ‘passive  bearing’  represented  a  plate 
clamped  and  loaded  at  its  two  remote  ends.  The  plate  had  a  hole  in  the  middle  which  contained  a 
non-load-bearing  bolt.  The  other  case  of  load  transfer,  referred  to  as  ‘active  bearing’,  represented 
a  plate  clamped  at  one  remote  end  and  loaded  through  a  pin  which  passed  through  a  hole  close  to 
the  non-clamped  end. 

The  base  laminate  consisted  of  an  alternating  0/90°  stacking  sequence  of  thirty  plies  made 
of  0.025"  thick  WR  glass  fabric  and  510-A  vinyl  ester  resin.  This  lay-up  was  fabricated  by 
Hardcore  DuPont  Composites,  L.L.C.,  into  several  panels  measuring  2  feet  by  3  feet  using  the 
SCRIMP  process.  The  post  cure  cycle  for  the  panels  was  2  hours  at  160°  F.  The  panels  were  then 
cut  into  specimens.  Figures  (la)  and  (lb)  show  the  geometry  of  typical  passive-bearing  and  active¬ 
bearing  specimens. 

To  investigate  the  behavior  of  composite  joints  under  static,  constant  amplitude,  and 
variable  amplitude  cyclic  loadings,  three  different  joint  configurations  were  designed  and 
fabricated.  These  three  configurations  were:  bonded  only,  bolted  only,  and  both  bonded  and 
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bolted.  Each  configuration  represented  a  single  butt-strap  joint.  Specimens  were  made  of  the 
laminate  described  above,  but  with  an  additional  fifteen  ply  laminate,  of  the  same  lay-up,  attached 
to  either  side  of  the  thicker  laminate.  Both  thicker  base  laminate  and  the  thinner  joint  adherends 
were  oriented  and  loaded  along  the  0/90°  material  direction.  Figure  (2)  shows  the  typical 
geometries  of  joint  specimens. 

STATIC  TESTS 

The  static  tests  consisted  of  ultimate  tension  and  compression  strength  tests  which  were 
designed  to  characterize  the  static  behavior  of  the  specimens  and  to  aid  in  selecting  the 
appropriate  stress  levels  for  the  fatigue  tests.  In  all  configurations,  some  test  specimens  were 
instrumented  with  strain  gages  in  order  to  monitor  the  structural  response  at  selected  points  of 
interest.  Strain  gage  locations  are  shown  in  Figures  (3)  through  (5).  Static  test  results  for  both 
specimens  and  joints  are  presented  in  Table  (1). 

The  static  test  results  for  the  passive  bearing  test  of  the  base  laminates  are  presented  in 
Figures  (6)  through  (14).  It  can  be  seen  fi'om  Figure  (6),  that  the  fracture  surfaces  of  tensile 
specimens  pass  through  the  bolt  hole  forming  a  90°  and  45°  angle  with  the  loading  direction,  for 
0/90°  and  ±45°  lay-ups  respectively.  For  the  compressive  tests,  the  ±45°  specimens  failed, 
forming  an  X-  shaped  damage  zone  centered  at  the  hole  and  swelling  in  the  intersection  of  the  two 
branches  of  the  X.  The  0/90°  compressive  specimen  was  loaded  up  to  the  capacity  of  the  testing 
machine  and  did  not  fail.  During  each  test,  applied  load  and  actuator  movement  were  recorded. 
Actuator  movement  was  measured  using  a  linearly  varying  displacement  transducer  (LVDT).  The 
load-deflection  curves  of  the  0/90°  specimens  show  a  significant  qualitative  difference  with  those  of 
the  ±45°  lay-up  (compare  Figure  (7)  with  Figure  (11),  and  Figure  (9)  with  Figure  (13)).  In  the 


former  figures,  the  curves  remain  almost  linear  up  to  the  failure  point,  while  in  the  latter  figures, 
there  is  a  considerable  increase  in  the  compliance  as  the  load  increases.  Moreover,  the  final 
deflection  reading  is  much  higher  than  that  of  the  cross  ply  counterpart.  This  is  obvious  because  of 
the  fundamental  difference  between  the  fiber-dominant  and  matrix-dominant  behavior  shown  by  the 
cross-ply  and  angle-ply  lay-ups,  respectively.  The  same  trend  can  be  observed  in  load-strain  curves 
(compare  Figure  (8)  with  Figure  (12),  and  Figure  (10)  with  Figure  (14)).  As  can  be  seen  on  the 
load-strain  curves,  the  constraint  provided  by  the  bolts  causes  the  compressive  strain  at  position 
number  4  of  the  passive  bearing  case  to  almost  vanish. 

The  static  test  results  for  the  active  bearing  case  are  shown  in  Figures  (15)  through  (25).  As 
shown  in  Figure  (15a),  in  tension  tests,  both  lay-ups  exhibited  a  shear-out  failure;  however,  the 
shear-out  surfaces  in  the  angle  ply  case  are  parallel  to  the  loading  direction,  while  those  of  the 
cross-ply  lay-up  are  inclined.  This  means  that  both  specimens  failed  in  the  least  fiber-dominant 
direction.  Unlike  the  passive  bearing  specimens,  the  stress-deflection  curves  for  the  tensile  test  of 
the  cross-ply  and  angle-ply  lay-ups  have  barely  any  qualitative  difference,  except  for  some  more 
ductility  in  the  0/90°  case,  as  shown  in  Figure  (16).  Regarding  strains,  the  somewhat  inconsistent 
final  portions  of  the  strain  versus  nominal  stress  traces,  may  be  due  to  local  differences  intensified 
by  the  localized  damage.  In  the  compression  tests,  the  stress-deflection  curves  are  totally  different. 
While  the  angle-ply  laminates  respond  almost  monotonically  to  failure,  the  cross-ply  laminate 
specimens,  besides  their  higher  strength,  show  a  step-like  trend,  meaning  that  some  kind  of 
secondary  resistance  forms  after  each  incremental  shear-out.  This  secondary  resistance  may  have 
been  produced  by  the  accumulation  of  the  material  removed  during  the  shear-out  process  in  fi-ont 
of  the  shear-out  path,  which  increases  the  effective  thickness.  Both  for  the  cross-ply  and  angle-ply 
lay-ups,  the  compressive  tests  were  terminated  due  to  the  limitations  dictated  by  the  fixture  and 
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specimen  configurations.  It  is  obvious  from  Figure  (15)  that  the  failure  is  only  a  partial  shear-out, 
indicating  that  both  specimens  may  have  carried  a  higher  load  if  the  testing  equipment  were  of 
greater  capacity.  However,  it  should  also  be  kept  in  mind  that  the  active  bearing  test  was  intended 
to  characterize  the  static  and  fatigue  bearing,  not  shear-out,  failure;  meaning  that  in  both  the 
tensile  and  compressive  strength  cases,  the  tests  could  be  stopped  before  actual  shear  out  occurred. 
In  fact,  in  all  static  active  bearing  case  tests,  bearing  failure  occurred  well  before  shear-out  failure. 
In  such  cases,  a  deviation  from  linearity  in  the  stress-deflection  curves  was  evidenced. 

The  static  test  results  of  the  joints  are  presented  in  Figures  (26)  through  (37).  Although  no 
photograph  is  available,  the  adhesive  failed  in  the  bonded  specimen  tests.  Figure  (26)  shows  the 
bolted-bonded  joint;  after  first  demonstrating  a  bond  failure,  the  specimen  finally  fails  in  a  shear  out 
mode.  Figure  (27)  shows  that  the  bolted  joint  exhibits  a  combination  of  cleavage-tension  and 
shear-out  failure.  In  the  compressive  regime  (Figures  (28)  and  (29)),  the  bonded  joint  has 
delaminated  in  both  laps  and  one  of  the  adherends  failed  due  to  the  combined  effect  of  contraction 
and  buckling.  Similarly,  the  bolted  joint  has  shortened  and  buckled  asymmetrically  (mode  2)  and 
all  joint  components  have  delaminated  as  a  consequence.  The  bonded-bolted  joint  showed  an  early 
bond  failure  and  a  final  failure  due  to  a  symmetric  buckling  and  a  consequent  delamination  in  all 
components.  Comparing  the  stress-deflection  curves  of  different  joints,  one  can  conclude  that  the 
bolted-bonded  joint  is  the  most  advantageous  configuration  as  far  as  the  tensile  strength  is 
concerned.  The  bolted-bonded  joint  had  a  higher  initial  stiflftiess,  a  higher  failure  load  and  a  very 
good  ductility  compared  with  the  other  two  specimens.  On  the  other  hand,  in  the  compressive 
regime,  although  it  had  a  higher  initial  stififtiess,  its  strength  and  ductility  were  lower  than  those  of 
the  bolted  configuration.  This  may  be  due  to  the  difference  in  the  buckling  mode,  since  a  mode  II 
buckling  requires  a  higher  load  in  order  to  occur. 
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FATIGUE  TESTS 


The  fatigue  tests  typically  consisted  of  two  or  three  specimens  tested  at  each  desired 
constant  amplitude  or  variable  amplitude  (RMS)  stress  level.  All  tests  were  performed  at  constant 
frequency  (the  frequency  and  other  pertinent  test  data  and  results  for  all  specimens  are  listed  in 
Tables  (2)  through  (7)).  The  variable  amplitude  loading  represents  a  segment  of  a  sinusoidal  time 
history  composed  of  10,000  end  points,  which,  with  the  assumption  of  stationarity  and  ergodicity, 
can  be  shown  to  be  a  narrowband  Gaussian  process.  This  segment  is  repeated  until  failure  occurs. 

The  data,  presented  here  for  each  specimen  type,  consist  of  experimental  fatigue  lives,  a 
brief  explanation  of  the  failure  modes,  peak  deflection  data,  temperature  data,  S-N  curves,  and  the 
random  fatigue  lives  predicted  by  the  Rayleigh  approximation.  A  comparison  between  the 
geometric  mean  of  the  experimental  fatigue  lives  and  predicted  fatigue  lives  can  be  found  in  Table 
(8).  Before  discussing  the  data,  the  significance  of  deflection  and  temperature  data  is  addressed. 

The  deflection  (actuator  displacement)  data  were  used  as  a  measure  of  the  overall  specimen 
compliance.  These  data  were  also  used  as  an  indicator  of  fatigue  damage.  Therefore,  when 
experimentally  investigating  the  damage  accumulation,  it  is  useful  to  collect  deflection  data.  Based 
on  the  peak  deflection  data,  the  evolution  of  damage  during  the  fatigue  life  can  be  calculated  from: 

D„=l-5i/5„  (47) 

where  Dn ,  5i,  and  6n  stand  for  the  deflection-based  damage  at  cycle  n,  initial  peak  deflection,  and 
the  peak  deflection  at  cycle  n,  respectively.  This  methodology  may  be  able  to  predict  the  variable 
amplitude  damage  accumulation  as  well. 
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The  acquisition  of  temperature  data  was  prompted  both  by  the  early  observations  in  the 
current  fatigue  test  program  and  by  the  survey  of  existing  literature.  During  the  fatigue  tests,  the 
temperature  of  the  specimens  increased  significantly  from  the  ambient  temperature  in  the  testing 
laboratory.  The  response  of  the  material  during  the  temperature  rise  was  so  severe,  especially 
during  the  low  cycle  fatigue  tests  with  relatively  high  frequency  (e.g.,  3  HZ),  that  the  material 
appeared  to  smolder,  producing  ash  and  a  burnt  plastic  smell  as  the  test  went  on.  Since  such  a 
reaction  is  clear  indication  of  a  change  in  mechanical  properties,  one  could  easily  conclude  that  the 
material  behavior  has  some  level  of  temperature  dependency.  This  conclusion,  when  accompanied 
with  the  creep  observed  during  the  static  tests,  indicated  that  the  material  was  visco¬ 
thermoelastic.  A  literature  search  revealed  that  the  temperature  dependence  of  FRPs  is  a  common 
concern  among  researchers.  It  has  been  observed  that  the  temperature  change  during  non- 
isothermal  conditions,  caused  by  thermal  or  mechanical  sources,  alters  the  material  behavior  and 
thereby  controls  the  mechanical  characteristics.  The  temperature  rise  may  be  caused  by  different 
damage  and/or  damping  processes  occurring  in  the  plastic  constituents,  especially  under  cyclic 
loads.  The  cyclic  energy  induced  into  the  system  is  dissipated  by  heat  convection,  different  forms 
of  damage  in  the  material,  and  hysteretic  damping.  The  latter  two  forms  of  energy  dissipation  are 
associated  with  temperature  increase  which,  in  turn,  deteriorates  the  mechanical  properties. 
Therefore,  the  mechanical  temperature  increase  and  its  consequences  are  expected  to  worsen  when 
the  load  frequency  is  increased.  Many  authors  have  investigated  the  temperature  effects  in  FRPs 
[Miyano  etal  (1994)  and  Nakai  and  Yamashita  (1994)].  Since  the  temperature  dependency  of  these 
materials  is  largely  due  to  the  properties  of  the  resin  matrix,  which  typically  has  a  higher  volume 
content  in  marine  type  FRPs,  the  consequences  of  temperature  dependence  can  be  more  severe  for 
marine  composites.  As  discussed  earlier,  the  plastic  constituents  of  FRPs  are  visco-thermoelastic 
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materials.  Therefore,  temperature  effects  are  coupled  with  creep,  relaxation,  and  time-under-load 
effects.  More  precisely,  while  the  temperature  effects  dominate  under  high  frequencies,  creep 
effects  dominate  at  low  frequencies.  In  their  1979  work.  Sun  and  Chan  [1979]  investigated  the 
interaction  of  temperature  and  creep  under  different  stress  levels  for  graphite  epoxy  composites. 
They  concluded  that  the  fatigue  strength  of  FRPs  varies  in  a  two  dimensional  space  formed  by 
stress  level  and  load  frequency.  The  latter  axis  reflects  the  effects  of  creep  and  temperature  in  low 
and  high  frequency  ranges,  respectively.  While  the  temperature  dependence  is  discussed  by  many 
authors  (e.g.,  Oberbach  and  Heese  [1972],  Nakada  et  al  [1991]),  many  others  have  considered  the 
other  extreme  of  the  frequency  axis,  investigating  the  time  and  creep  dependence  of  the  fatigue 
properties  of  FRPs.  Mandell  and  Meier  [1983]  have  investigated  the  effect  of  loading  time  on  the 
tensile  fatigue  of  glass-reinforced  epoxy.  Morinaka,  Fujii,  and  Amijima  [1991]  have  studied  the 
effect  of  loading  and  unloading  duration  on  the  fatigue  damage  process  of  plain  woven  glass  fiber 
reinforced  plastics.  They  concluded  that  the  fatigue  life  is  affected  by  both  the  loading  and 
unloading  duration.  Based  on  these  and  similar  observations,  which  is  not  our  goal  to 
comprehensively  cover,  it  seems  necessary  to  incorporate  the  temperature,  frequency  and  creep 
effects  in  fatigue  life  and  damage  accumulation  models.  In  an  effort  to  incorporate  the  temperature 
and  creep  effects  in  fatigue  life  models.  Sun  and  Chan  in  the  aforementioned  paper,  have  presented 
two  separate  fatigue  life  models  for  temperature-dominant  and  creep-dominant  behaviors. 
Amijima,  et  al.  [1982]  have  found  an  experimental  frequency  dependent  S-N  curve.  Tsai,  et  al 
[1987]  assumed  the  strain  to  beas  the  governing  variable  and  defined  the  dynamic  modulus  (which 
is  affected  by  test  frequency)  so  as  to  incorporate  the  frequency  effects  in  their  fatigue  life  model. 
Miyano  et  al  [1994]  have  implemented  the  invariability  of  the  slope  of  the  S-N  curves  with 
temperature  and  frequency  to  derive  a  relationship  for  predicting  the  flexural  fatigue  strength  at 
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arbitrary  conditions.  While  the  aforementioned  fatigue  life  models  have  taken  into  account  the 
temperature  dependence  of  damage,  it  is  important  to  note  that  the  material  temperature  is,  in  turn, 
damage  dependent.  In  particular,  the  damage  itself  can  cause  a  substantial  material  temperature 
increase.  Therefore,  logically,  temperature  can  be  used  as  a  measure  of  damage.  For  example, 
Neubert,  Schulte,  and  Harig  [1990]  have  investigated  the  damage  development  in  CFRPs  by 
monitoring  load-induced  material  temperature  changes.  However,  in  applying  the  change  in 
temperature  as  a  measure  of  damage,  one  should  be  extremely  cautious,  since  the  temperature  rise 
cannot  be  attributed  solely  to  the  damage.  Other  mechanisms,  such  as  hysteretic  damping,  may 
also  contribute  to  the  temperature  increase.  Therefore,  the  accurate  temperature-based  damage 
accumulation  model  is  the  one  in  which  other  effects  are  excluded.  The  following  sections  describe 
the  data  collected  for  each  of  the  passive,  active,  and  joint  specimens. 

PASSIVE  BEARING  CASE 

Experimental  Fatigue  Lives  ;  The  fatigue  lives  obtained  from  constant  and  variable 
amplitude  tests  of  the  passive  bearing  case  are  presented  in  Tables  (2)  and  (3),  respectively. 

Failure  Modes  ;  The  failure  for  all  fatigue  specimens  was  defined  as  occurring  when  the 
peak  deflection  reached  four  times  the  initial  value.  For  the  variable  amplitude  case,  this  definition 
considers  the  largest  of  the  peaks  in  each  of  10,000  load  points.  For  0/90"  lay-ups,  at  failure,  a 
transverse  whitened  damage  band  appeared  in  the  specimen  passing  across  the  bolt  hole,  as  shown 
in  Figure  (38b).  For  the  ±  45"  angle-ply  lay-up,  however,  damage  was  in  the  form  of  an  X-shaped 
whitened  band  passing  through  the  gage  section,  with  each  branch  of  X  forming  approximately  a 
45“  angle  with  the  loading  direction,  as  shwon  in  Figure  (38a). 

Deflection  Data  :  The  deflection  data  are  presented  in  the  form  of  peak  deflection  versus 
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number  of  cycles  curves  in  Figures  (39)  through  (41).  All  these  figures,  and  subsequent  figures, 
indicate  that  the  tensile  deflection  governs  the  failure;  and,  in  most  cases,  the  compressive 
deflection  does  not  change  considerably.  The  obvious  trend  of  the  tensile  deflection  curves 
indicates  that,  no  matter  what  the  stress  level  or  lay-up,  the  peak  deflection  does  not  change 
considerably  until  the  final  failure  is  approached.  According  to  the  discussion  in  the  previous 
paragraphs,  a  compliance-based  experimental  damage  accumulation  curve  can  be  obtained  fi’om 
equation  (47).  The  results  for  some  specimens  of  the  cross-ply  and  angle-ply  lay-ups  are  shown  in 
Figures  (42)  and  (43),  respectively. 

Temperature  Data  :  The  temperature  data  are  shown  in  the  form  of  temperature  versus 
normalized  number  of  cycles  in  Figures  (44)  through  (52)  for  different  stress  levels,  lay-ups,  and 
loading  scenarios.  The  large  number  of  specimens  shows  the  wide  qualitative  difference  between 
the  curves  of  different  cases,  and  indicates  the  temperature  rise  is  a  consequence  of  interaction  of 
many  parameters.  The  specimen  temperature  can  build  up  in  different  ways,  depending  on  how 
those  parameters  interact.  Another  noteworthy  point  is  that  there  is  a  close  resemblance  between 
Figures  (42)  and  (48)  and  between  Figures  (43)  and  (47);  which  means  that  the  temperature  rise 
has  been  able  to  capture  the  main  features  of  the  compliance-based  damage  accumulation  curves 
very  well.  This  may  indicate  that  temperature  is  an  alternative  basis  for  phenomenological  damage 
accumulation  models. 

S-N  Curves  :  The  experimental  fatigue  lives  from  the  constant-amplitude  tests  of  the 
passive-bearing  specimens  are  presented  in  Figure  (53)  on  a  log(S)-log(N)  scale  and  in  Figure  (54) 
on  a  S-log(N)  scale.  Both  logarithmic  and  semi-logarithmic  S-N  curve  models  were  fit  to  these 
data  by  linear  regression.  It  can  be  seen  that  the  experimental  data  from  the  0/90°  specimens 
exhibit  very  little  scatter,  and  although  the  data  fit  both  S-N  curve  models  fairly  well,  a  slightly 
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better  fit  can  be  seen  using  the  Log-Log  model.  The  data  from  +  45°  test  setup  exhibit  much  more 
scatter,  and  although  not  fitting  either  S/N  curve  model  particularly  well,  appear  to  fit  the 
logarithmic  model  better  than  the  semi-logarithmic  model. 

Random  Amplitude  Fatigue  Lives  :  Experimental  variable  amplitude  fatigue  lives  are 
obtained  by  applying  the  aforementioned  time  history  to  the  specimens  at  different  RMS  stress 
levels;  and  the  analytic  fatigue  lives  are  calculated  by  following  the  fatigue  damage  accumulation 
models  (Equations  33  and  38)  discussed  previously.  Comparison  between  the  variable-amplitude 
fatigue  results  and  the  Rayleigh  approximation,  based  on  both  logarithmic  and  semi-logarithmic  S- 
N  curve  models,  is  provided  in  Table  (8).  It  can  be  seen  that,  for  the  ±45°  specimens  and  for 
both  RMS  stress  levels  considered,  the  classical  Rayleigh  approximation  underestimates  the  actual 
fatigue  life  by  about  60%.  The  same  trend  is  observed  for  0/90°  specimens  for  the  higher  RMS 
stress  level  of  6.0  ksi.  At  the  same  time,  the  fatigue  life  at  the  lower  RMS  stress  level  of  4.5  ksi  is 
significantly  overestimated  by  Equation  (33).  For  both  specimen  types,  the  predictions  at  the 
higher  RMS  stress  levels  (2.5  and  6.0  ksi,  respectively)  are  somewhat  improved  by  the  use  of 
Equation  (38),  but  this  improvement  is  quite  marginal.  At  the  same  time.  Equation  (38)  leads  to 
further  deviation  from  the  experimental  mean  for  the  two  lower  RMS  stress  levels  (1.75  and  4.5 
ksi,  respectively).  These  large  deviations  mainly  come  from  the  basic  assumptions  built  in  the 
Rayleigh  approximation. 

As  pointed  out  in  the  previous  section,  the  classical  Rayleigh  approximation  formula. 
Equation  (33),  is  based  on  both  Equation  (1)  and  on  the  Palmgren-Miner  damage  accumulation 
rule  Equation  (29).  Therefore,  misrepresentation  of  the  actual  stress-life  relationship  caused  by  a 
poor  S-N  curve  function,  insufficient  S-N  data,  or  disregarding  the  sequence  effects  by  using  the 
Miner  rule,  are  main  reasons  for  the  shortcomings  of  the  Rayleigh  approximation  in  estimating  the 
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variable  amplitude  fatigue  lives. 


ACTIVE  BEARING  CASE 

Experimental  Fatigue  Lives  :  The  constant  amplitude  and  variable  amplitude  fatigue  lives 
for  the  active  bearing  case  are  presented  in  Tables  (4)  and  (5),  respectively. 

Failure  Modes  :  The  active  bearing  specimens  were  intended  to  characterize  the  bearing 
fatigue  behavior  and  failure  of  the  composite  material  under  constant  and  variable  amplitude 
loading.  The  failure  was  defined  as  occurring  when  the  peak  deflection  reached  four  times  the 
original  peak  deflection.  Figure  (55)  shows  the  typical  appearance  of  cross-ply  and  angle-ply 
specimens  after  failure.  Obviously,  there  is  no  difference  between  the  bearing  failure  modes  of  the 
cross-ply  and  angle-ply  specimens.  The  elliptic  glassy  region  circumscribing  the  hole  for  both  lay¬ 
ups  seems  to  be  a  consequence  of  local  delamination  which  occurs  due  to  bearing  stresses. 

Peak  deflection  data  ;  Figures  (56)  through  (58)  show  the  evolution  of  the  peak  deflection 
during  the  constant  amplitude  fatigue  tests,  for  the  active  bearing  case.  Unlike  the  passive  case,  the 
curves  have  a  shorter  low  cycle  flat  portion,  and  then  suddenly  the  deflection  begins  to  build-up 
roughly  midway  through  the  fatigue  life.  Also,  the  downward  curvature  of  the  final  part  of  the 
curves  distinguish  the  active  bearing  curves  from  those  of  the  passive  case.  A  common  feature  of 
all  cases,  however,  is  that  the  failure  is  controlled  by  the  tensile  deflection. 

Temperature  data  :  Figures  (59)  through  (66)  present  the  temperature  rise  during  the 
fatigue  tests  of  active  bearing  case.  These  plots  do  not  exhibit  the  same  behavior  as  the  passive 
bearing  plots  nor  are  they  as  diverse.  Some  of  the  bearing  case  plots  show  a  short  initial  increase 
and  some  do  not.  Prior  to  the  final  build  up,  some  of  the  active  bearing  curves  show  almost  no 
variation  in  the  temperature,  while  for  some  others  there  is  some  temperature  variation.  Also,  the 
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duration  of  the  final  temperature  rise  differs  considerably  for  different  specimens. 

S-N  Curves  :  The  S-N  curves  for  the  active  bearing  specimens  are  shown  in  Figures  (67) 
and  (68)  in  log-log  and  semi-log  scales,  respectively.  From  the  figures,  it  can  be  seen  that  the 
scatter  of  the  cross-ply  data  is  less  than  that  of  the  angle  ply  data.  This  agrees  with  the  passive 
bearing  case  data. 

Random  Fatigue  Lives  :  Table  (8)  summarizes  the  random  fatigue  lives  obtained  from  the 
experiments,  as  well  as  the  calculated  lives.  Generally,  the  difference  between  the  two  lives  is  larger 
for  lower  stress  levels.  Like  the  passive  bearing  case,  the  large  difference  between  the 
experimental  and  theoretical  results  is  not  consistently  conservative,  or  non-conservative.  Effects 
such  as  the  quality  with  which  the  constant  amplitude  data  are  represented,  load  sequence  and 
load  frequency  may  very  well  be  contributing  to  this  lack  of  consistency. 

JOINTS 

Experimental  Fatigue  Lives  :  Tables  (6)  and  (7)  show  the  constant  amplitude  and  variable 
amplitude  test  results  for  different  types  of  composite  joints,  respectively. 

Failure  Modes  :  Figures  (69)  and  (70)  show  the  typical  failures  of  bonded  only,  bolted  only, 
and  bolted-bonded  joints.  The  failure  of  the  bonded  joints  is  an  adhesive  failure  of  the  bond,  while 
that  of  the  bolted  joints  is  a  bearing  followed  by  a  partial  shear-out  failure.  The  white  powder 
around  the  bolt  holes  is  a  result  of  the  frictional  interaction  of  the  parts.  The  failure  of  the  bonded- 
bolted  joints  occurs  in  two  stages.  First,  and  very  early,  comes  the  adhesive  failure  of  the  bonds 
and  then,  typically  after  a  long  time,  the  bearing  and  partial  shear-out  of  the  bolts  occurs. 

Deflection  Data  :  Figures  (71)  through  (73)  show  the  peak  deflection  versus  number  of 
cycles  curves  for  the  bonded  joints.  Note  that  the  curves  in  Figures  (72)  and  (73)  are  enlarged  in 
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the  final  failure  part  to  help  distinguish  the  step-like  growth  of  deflection  to  failure.  However,  such 
a  trend  was  not  observed  in  Figure  (71).  This  step-wise  failure  is  not  of  significance  when  studying 
the  fatigue  damage  accumulation  and  life.  However,  it  may  be  helpful  in  failure  mechanism 
investigations.  Figures  (74)  through  (77)  show  the  evolution  of  deflection  for  bolted  joints,  and 
Figures  (78)  through  (80)  for  bolted-bonded  joints.  From  Figures  (74)  through  (77),  the  trend 
resembles  that  of  the  passive  bearing  case,  except  that  in  Figure  (77),  the  final  build  up  is  more 
gradual.  In  the  deflection  curves  of  the  bolted-bonded  joints,  the  final  sudden  build  up  seems  to 
follow  a  small  step-like  growth.  Figures  (81)  through  (83)  exemplify  the  deflection-based  damage 
accumulation  for  different  types  of  joints.  It  is  interesting  to  note  the  fundamental  difference 
between  the  three  types  of  joints.  For  bonded  joints,  the  stress  level  has  a  marginal  effect  on  the 
damage  accumulation  form.  On  the  other  hand,  the  damage  accumulation  in  bolted  joints  is 
governed  by  the  stress  level.  The  bolted-bonded  specimens  have  shown  an  intermediate  trend. 

Temperature  Data  :  Figure  (84)  shows  an  example  of  temperature  change  in  a  bolted- 
bonded  composite  joint.  As  expected,  there  is  a  considerable  difference  between  the  temperature  of 
the  washer  and  that  of  the  composite  material  in  the  vicinity  of  the  washer. 

S-N  Curves  :  The  constant-amplitude  data  for  all  joint  configurations  are  summarized  in  the 
form  of  S-N  curves  in  Figure  (85)  and  (86).  The  curves  indicate  a  definite  deviation  from  the 
logarithmic  S-N  curve  model  for  all  joint  configurations.  Neither  S/N  curve  model  appears  to 
represent  the  data  particularly  well. 

Random  Lives  :  The  experimental  and  theoretical  random  lives  are  summarized  in  Table 
(8).  Like  the  previous  two  cases,  it  can  be  seen  that  the  Rayleigh  approximation  generally  has  not 
been  able  to  predict  the  lives  accurately.  Predictions  tend  to  be  conservative  low  RMS  levels  and 
non-conservative  at  higher  RMS  levels.  This  is  most  likely  due  to  the  poor  representation  of  the 
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S/N  curves  to  constant  amplitude  data,  or  possibly  due  to  sequence  effects. 


CONCLUSIONS 


Fatigue  damage  in  FRP  laminates  occurs  in  the  form  of  diffuse  zones  containing  micro¬ 
cracks  in  the  matrix,  fiber-matrix  debonding  sites,  ply  delamination,  and  broken  fibers.  This  damage 
leads  to  degradation  of  the  elastic  moduli  of  the  laminate;  FRP  fatigue  damage  accumulation  is  best 
modeled  by  stiflftiess  degradation  approaches.  Various  S-N  curve  models  for  constant-amplitude 
fatigue  damage  characterization  can  be  described  in  terms  of  the  fatigue  modulus  concept  of 
Hwang  and  Han.  Very  little  investigation  has  been  done  using  fatigue  damage  accumulation 
techniques.  Most  available  data  deal  with  two-stage  loading  cases  only.  A  few  cumulative  damage 
hypotheses  have  been  proposed.  However,  due  to  limited  experimental  data,  it  is  not  clear  which 
of  these  hypotheses  are  most  suitable  for  naval  applications.  The  constant  and  variable  amplitude 
fatigue  studies  on  WR  laminates  and  joints  conducted  at  The  George  Washington  University 
indicate  a  fair  fit  of  the  S-N  data  for  the  cross-ply  laminate  of  both  passive  and  active  bearing  cases 
to  the  power  law  S-N  curve.  The  fit  is  not  as  good  for  the  angle-ply  lay-ups.  At  the  same  time, 
the  laminate  joints  tested  show  a  definite  trend  toward  deviation  from  the  power  law  S-N  curve 
model.  Random-amplitude  fatigue  tests  with  narrowband  stress  histories  demonstrate  significant 
deviation  of  the  observed  fatigue  lives  from  those  predicted  with  the  Rayleigh  approximation; 
tending  to  overpredict  for  joints  at  the  high  RMS  level  and  underpredict  for  the  low  RMS  level. 
Although  not  as  consistent,  the  active  and  passive  bearing  specimens  tend  to  exhibit  just  the 
opposite  behavior,  overpredicting  at  the  low  RMS  level  and  underpredicting  at  the  high  RMS  level. 
Such  agreement  may  be  attributed  to  an  inadequacy  of  the  Palmgren-Miner’s  damage  accumulation 
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rule  to  accurately  reflect  the  fatigue  damage  accumulation  behavior  in  composite  materials.  Other 
reasons  for  the  deviation  include  a  poor  representation  or  insufficiency  of  the  constant  amplitude 
data  to  be  represented  by  an  analytical  expression,  as  well  as  temperature  and  load  frequency 
effects  which  have  not  yet  been  quantified  and  incorporated  into  a  fatigue  damage  accumulation 
model.  In  general  however,  use  of  the  Log-Log  S/N  curves  produced  slightly  better  predictions. 

RECOMMENDATIONS 


The  data  generated  and  presented  in  this  report  provide  an  initial  step  in  understanding 
fatigue  damage  accumulation  in  marine-grade  composite  laminates  and  joints  subjected  to  constant 
and  random  amplitude  loadings.  Several  problems  have  been  identified  in  the  course  of  this  effort; 
most  notably,  the  inability  to  properly  represent  the  constant  amplitude  fatigue  data  by  traditional 
power  law  forms  and  a  general  lack  of  sufficient  data.  This  problem  alone  may  have  been  the  most 
important  factor  in  the  inability  to  accurately  and  consistently  predict  failure  lives  of  test  specimens 
subjected  to  random  loads.  Other  effects,  such  as  temperature  and  load  frequency,  undoubtedly 
have  influenced  these  results,  but  it  is  not  known  by  how  much. 

Follow-on  efforts  should  include  characterizing  temperature  and  load  frequency  quantities 
so  their  effects  can  be  incorporated  into  a  fatigue  damage  accumulation  model.  In  addition, 
predictions  should  be  made  with  the  other  models  described  in  this  report  and  compared  with  the 
experimental  results.  Further,  tests  under  constant  and  random  amplitude  loadings  will  likely  need 
to  be  repeated,  this  time  varying  load  frequency  while  monitoring  the  temperature. 

Different  geometries  and  methods  of  attachment/adhesives  should  also  be  investigated  to 
improve  the  fatigue  behavior  of  both  laminates  and  joints. 
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3  3/4 


(a) 

Typical  Bonded  only  Specimen 


Figure  2  -  Geometry  of  Different  Types  of  Composite  Joints 
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Gage 
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Figure  3  -  Strain  Gage  Locations  in  Passive  Bearing  Case;  Static  Test  Specimens 
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Figure  4  -  Strain  Gage  Locations  in  Active  Bearing  Case;  Static  Test  Specimens 


17.75" 


Strain  Gage  Positions  for  Bolted  only  and  Bolted  &  Bonded  Specimens 


17.75" 


Strain  Gage  Positions  for  Bonded  only  Specimens 


Figure  5  -  Strain  Gage  Locations  in  Static  Tests  of  Joint  Specimens 
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Figure  6  -  Failure  Modes  for  Strength  Tests  of  Passive  Bearing  Case 


ecimen  6-0750-6  (Ultimate  Tensile  Strength  =  37,500  psi) 
ecimen  6-0750-7  (Ultimate  Tensile  Strength  =  38,000  psi) 
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Figure  7  -  Stress-Deflection  Curves  for  Ultimate  Tensile  Strength  Tests  of  0/90°  Passive  Specimens 


42000 
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Figure  8  -  Strain  Gage  Readings  During  Ultimate  Tensile  Strength  Tests  of  0/90"  Passive  Specimens;  Specimen  6-0750-6 
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Figure  9  -  Stress-Deflection  Curve  for  Ultimate  Compressive  Strength  Tests  of  0/90"  Passive  Specimen 
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Figure  10  -  Strain  Gage  Readings  During  Compressive  Strength  Test  of  0/90“  Passive  Specimens 
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Figure  1 1  -  Stress-Deflection  Curves  for  Ultimate  Tensile  Strength  Tests  of  ±  45°  Passive  Specimens 
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Gage  Readings  During  Ultimate  Tensile  Strength  Tests  of  ±45”  Passive  Specimens;  Specimen  9-0750-9 
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■Deflection  Curves  for  Ultimate  Compressive  Strength  Tests  of  ±45"  Passive  Specimens 
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Figure  14  -  Strain  Gage  Readings  During  Ultimate  Compressive  Strength  Tests  of  ±45°  Passive  Specimens;  Specimen  9-0750-8 
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Figure  15  -  Failure  Modes  for  Strength  Tests  of  Active  Bearing  Case 


13000 
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Figure  16  -  Stress-Deflection  Curves  for  Active  Bearing  Case:  Tensile  Strength  Tests 


62 


Gage  Readings  for  Tensile  Strength  Tests  of  Active  Bearing  Specimens;  Specimen  9-0750-4  (±  45“) 
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Gage  Readings  for  Tensile  Strength  Tests  of  Active  Bearing  Specimens;  Specimen  9-0750-5  (±45°) 
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Figure  19  -  Strain  Gage  Readings  for  Tensile  Strength  Tests  of  Active  Bearing  Specimens;  Specimen  6-0750-10  (0/90°) 
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Figure  20  -  Strain  Gage  Readings  for  Tensile  Strength  Tests  of  Active  Bearing  Specimens;  Specimen  6-0750-9  (0/90°) 
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Figure  21  -  Stress-Deflection  Curves  for  Compressive  Strength  Tests  of  Active  Bearing  Case 
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Figure  22  -  Strain  Gage  Readings  in  Active  Bearing  Case:  Compressive  Strength  Tests;  Specimen  9-075 


ressive  Strains 
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Figure  23  -  Strain  Gage  Readings  in  Active  Bearing  Case:  Compressive  Strength  Tests;  Specimen  9-075 
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Figure  24  -  Strain  Gage  Readings  in  Active  Bearing  Case:  Compressive  Strength  Tests;  Specimen  6-0750-12  (0/90°) 


pressive  Strains 
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Figure  30  -  Stress-Deflection  for  Tensile  Strength  Tests  of  Composite  Joints 
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Figure  31  -  Stress-Deflection  for  Compressive  Strength  Tests  of  Composite  Joints 
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Figure  35  -  Strain  Gage  Readings  for  Compressive  Strength  Test  of  Bolted  Joints 


00091 


81 


32000 


82 


84 


Figure  39  -  The  Evolution  of  Deflection  During  the  Fatigue  Life  of  Passive  Specimen  6-0750-1:  0/90",  7.5  ksi 
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Figure  40  -  The  Evolution  of  Deflection  During  the  Fatigue  Life  of  Passive  Specimen  6-0750-8;  0/90°,  15  ksi 
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Figure  41  -  The  Evolution  of  Deflection  During  the  Fatigue  Life  of  Passive  Specimen  9-0750-7:  ±  45“,  7.5  ksi 
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-  Constant  Amplitude  Deflection-Based  Fatigue  Damage  Accumulation  in  0/90“  FRP  Laminates 
with  Holes  -  Passive  Bearing  Case 
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-  Constant  Amplitude  Deflection-Based  Fatigue  Damage  Accumulation  in  ±  45°  FRP  Laminates 
with  Holes  -  Passive  Bearing  Case 


Frequency:  3  Hz 
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Temperature  Rise  During  Constant  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  Specimen  8-0750-3.  ±  45 , 3.5  ksi 


Failed  after  6,947,424  Cycles 
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Temperature  Rise  During  Constant  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  Specimen  9-0750-5;  ±45°,  3.5 
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Figure  46  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  Specimen  9-0750-6:  ±45°,  5  ksi 
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Temperature  Rise  During  Constant  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  ±  45°,  7.5  ksi 
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Figure  48  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  Specimen  6-0750-1;  0/90°,  7.5  ksi 


Specimen  6-0750-3  (Failed  after  474.871  Cycies) 
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Figure  49  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  0/90°,  10  ksi 
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Frequency:  1  Hz 
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Figure  51  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  0/90°,  20  ksi 


130 
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Figure  52  -  Temperature  Rise  During  Variable  Amplitude  Fatigue  Tests  of  Passive  Bearing  Case;  0/90°,  RMS=4.5  ksi 
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Figure  53  -  S-N  Curves  for  Passive  Bearing  Case  (Log-Log  Curve  Fit) 


0/90  Degree  Specimens 
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Figure  54  -  S-N  Curves  for  Passive  Bearing  Case  (Semi-Log  Curve  Fit) 
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Figure  55  -  Fatigue  Failure  Mode  of  Active  Bearing  Specimens  (Both  Lay-ups) 
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Figure  56  -  Evolution  of  Deflection  During  the  Fatigue  Life  of  Active  Bearing  Specimens; 
Specimen  8-0750-2:  ±  45“,  2.5  ksi,  Constant  Amplitude 


Evolution  of  Deflection  During  the  Fatigue  Life  of  Active  Bearing  Specimens; 
Specimen  8-0750-4;  ±45°,  2.5  ksi.  Constant  Amplitude 
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Figure  58  -  Evolution  of  Deflection  During  the  Fatigue  Life  of  Active  Bearing  Specimens; 
Specimen  8-0750-1  ;  ±45°,  2.5  ksi.  Constant  Amplitude 
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Figure  59  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Testing  for  +  45°  FRP  Laminates 
with  Holes  (2.5  ksi)  -  Active  Bearing  Case  -  8-0750-1 


130 


105 


Figure  60  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Testing  for  ±  45  "  FRP  Laminates 
with  Holes  (2.5  ksi)  -  Active  Bearing  Case  -  8-0750-2 
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Figure  61  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Testing  for  ±  45  °  FRP  Laminates 
with  Holes  (2.5  ksi)  -  Active  Bearing  Case  -  8-0750-4 
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Figure  62  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  Testing  for  ±  45"  FRP  Laminates 
with  Holes  (3  ksi)  -  Active  Bearing  Case  -  7-0750-12 


Specimen  7-0750-14  -  Failed  after  44,405  Cycles 
Specimen  7-0750-1 1  -  Failed  after  43,807  Cycles 
Specimen  8-0750-5  -  Failed  after  55,221  Cycles 
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Figure  63  -  Temperature  Rise  During  Variable  Amplitude  Fatigue  Testing  for  ±  45"  FRP  Laminates 
with  Holes  (1.75  ksi)  -  Active  Bearing  Case 
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Figure  64  -  Temperature  Rise  During  Variable  Amplitude  Fatigue  Testing  for  ±  45”  FRP  Laminates 
with  Holes  (1.25  ksi)  -  Active  Bearing  Case 
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Figure  65  -  Temperature  Rise  During  Variable  Amplitude  Fatigue  Testing  for  0/90°  FRP  Laminates 
with  Holes  (1 .75  ksi)  -  Active  Bearing  Case 


Ill 


Figure  66  -  Temperature  Rise  During  Variable  Amplitude  Fatigue  Testing  for  0/90”  FRP  Laminates 
with  Holes  (1.25  ksi)  -  Active  Bearing  Case 
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Figure  67  -  S-N  Curves  for  Active  Bearing  Case  (Log-Log  Curve  Fit) 


ecimens 

ecimens 

Fit-!og(N)  =  7.981  -  0.780Se 
Fit-Iog{N)  =  8.218  -  0.992  Se 
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Figure  68  -  S-N  Curves  of  Active  Bearing  Case  (Semi-Log  Curve  Fit) 


Figure  69  -  Appearance  of  Composite  Joints  after  Fatigue  Tests 
(FROM  LEFT  :  BONDED,  BOLTED,  AND  BONDED-BOLTED  JOINTS) 
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Figure  71  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bonded  Joints;  Specimen  JA-5:  5  ksi 


Tensile  LVDT  -  Failed  after  86,206  Cycles 
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Figure  72  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bonded  Joints;  Specimen  JA-18:  4  ksi 
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Figure  73  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bonded  Joints;  Specimen  JA-8:  3  ksi 
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Figure  74  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bolted  Joints;  Specimen  JB-5;  5  ksi 
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Figure  75  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bolted  Joints;  Specimen  JB-1 1:3.5  ksi 
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Figure  76  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bolted  Joints;  Specimen  JB-13:  3  ksi 
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Figure  77  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bolted  Joints;  Specimen  JB-14;  3  ksi 


Peak  Tensile  LVDT  -  Failed  after  1,481  Cycles 
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Figure  79  -  Evolution  of  Deflection  During  the  Constant  Amplitude  Fatigue  Tests  of  Bolted-Bonded  Joints;  Specimen  JAB-3;  7.5  ksi 
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iled  after  86,206  Cycles 
ed  after  3,167  Cycles 
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igure  81  -  Constant  Amplitude  Deflection-Based  Fatigue  Damage  Accumulation  in  Bonded  Joints 


Specimen  JB-3  -  Stress  Level  7.5  ksi  -  Failed  after  305  Cycles 

Specimen  JB-6  -  Stress  Level  5  ksi  -  Failed  after  799  Cycles 
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Figure  82  -  Constant  Amplitude  Deflection-Based  Fatigue  Damage  Accumulation  in  Bolted  Joints 
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Figure  83  -  Constant  Amplitude  Deflection-Based  Fatigue  Damage  Accumulation  in  Bolted-Bonded  Joints 
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Figure  84  -  Temperature  Rise  During  Constant  Amplitude  Fatigue  of  a  Bolted-Bonded  Joint;  Specimen  JAB-7:  6  ksi 


Bolted  &  Bonded 
Bonded  only 
Bolted  only 
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Figure  85  -  S-N  Curves  of  FRP  Laminate  Joints  (Log-Log  Curve  Fit) 


Bolted  &  Bonded 
Bonded  only 
Bolted  only 
Bolted  only  (Ran  out) 
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Figure  86  -  S-N  Curves  of  FRP  Laminate  Joints  (Semi-Log  Curve  Fit) 


TABLE  1  -  STATIC  TEST  RESULTS  FOR  COMPOSITE  SPECIMENS  AND  JOINTS 
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TABLE  2  -  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  FRP  LAMINATES  -  PASSIVE  BEARING  CASE 


Additional  Data 
Collected 

Temp. 

1 

1 

1 

1 

s 

S 

1 

1 

1 

s 

s 

s 

S 

1 

1 

1 

1 

LVDT 

1 

1 

1 
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1 

1 

S 
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1 

1 

1 

s 

s 

s 

1 

1 

1 

1 

Cycles  to 
Failure 

20,000,000  + 

19,100,000  + 

1,204,300 

6,946,300 

21,800 

o 

o 

vo" 

7,530 

490 

570 

540 

2,500 

240 

Failure  Criterion 
(x  Initial  LVDT) 

VD 

VO 

Load 

Frequency 

(Hz) 

* 

cn 

- 

ro 

in 

t-H 

Max. 

Load 

Level  (lbs) 

5,858 
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5,864 

8,112 

00 

oo'' 
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12,023 

12,424 

17,914 

17,924 

17,773 

24,000 
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Stress 

Level 

(ksi) 

2.5 

2.5 
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3.5 
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Identification 

No. 

9-0750-4 

8-0750-6 

8-0750-3 

9-0750-5 
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o 

o 

i 

00 

9-0750-1 

9-0750-6 

8-0750-4 

8-0750-11 

9-0750-7 

9-0750-2 

9-0750-3 

Lay¬ 

up 

+45° 
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+  Indicates  test  was  suspended  without  failure 

*Afler  1,969,795  cycles  the  load  frequency  was  changed  to  3  Hz. 

Note:  All  +45°  specimens  developed  X-Shaped  damage  bands  along  the  principal  material  directions  in  the  gage  section. 
These  damage  bands  may  be  due  to  shear  (mode  II)  cracking  of  the  laminates. 


TABLE  2  (con’t)  -  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  FRP  LAMINATES  -  PASSIVE  BEARING  CASE 
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TABLE  3  -  VARIABLE  AMPLITUDE  FATIGUE  TESTS  OF  FRP  LAMINATES  -  PASSIVE  BEARING  CASE 


Additional  Data 
Collected 

Temperature 

1 

1 

1 
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! 
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1 

1 

1 

1 

S 

S 

s 
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1 

S 

Cycles  to 
Failure 
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22,600 

857,500 

566,800 
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±7.597 

00 

d 

+1 

00 

d 

+I 

00 

d 

+I 

±19.536 

±19.536 

±19.536 

±26.048 

±26.048 

±26.048 

RMS 

(ksi) 

1.75 

1.75 

1.75 

2.5 

2.5 

2.5 

4.5 

4.5 

■ 

VO 

■ 

■ 

Identiflcation 

No. 

8-0750-2 

8-0750-8 

o 

T-H 

1 

o 

o 

1 

00 

1 

o 

o 

I 

00 

8-0750-5 

8-0750-9 

7-0750-2 

7-0750-3 

1 

o 

O 

1 

VO 

1 

O 

vn 

r- 

o 

1 

Lay-up 

±45° 

0/90° 
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TABLE  4  -  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  FRP  LAMINATES  -  ACTIVE  BEARING  CASE 


Additional  Data 

Temp. 

S 

S 

S 

S 

1 

1 

I 

1 

f 

1 

1 

1 

j 

LVDT 

S 

S 

B 

B 

B 

B 

B 

B 

B 

1 

1 

Cycles  to 
Failure 

238,000 

268,000 

340,900 

100,300 

344,700 

1,104,600 

24,600 

39,500 

- 1 

1,170 

970 

Failure  Criterion 
(x  Initial  LVDT) 

Load 

Frequency 

(Hz) 

- 

r— 1 

- 

1 

1 

1 

- 

- 

- 

Max.  Load 
Level  (lbs) 

5,790 

5,849 

5,867 

7,275 

1 _ 

7,103 

7,205 

9,910 

10,054 

12,323 

12,433 

Max.  Stress 
Level  (ksi) 

2.5 

2.5 

2.5 

m 

cn 

ro 

IT) 

in 

Identification 

No. 

8-0750-1 

8-0750-2 

8-0750-4 

7-0750-12 

8-0750-3 

9-0750-1 

9-0750-2 

o 

tT) 

r- 

o 

CN 

7-0750-9 

9-0750-3 

Lay-up 

±45° 
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TABLE  4  (con’t)  -  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  FRP  LAMINATES  -  ACTIVE  BEARING  CASE 
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TABLE  5  -  VARIABLE  AMPLITUDE  FATIGUE  TESTS  OF  FRP  LAMINATES  -  ACTIVE  BEARING  CASE 


Additional  Data 
Collected 

Temperature 

S 

s 

S 

S 

S 

S 

s 

s 

S 

s 

s 

Cycles  to 
Failure 

228,100 

198,100 

474,500 

43,800 

o 

o 

55,500 

231,700 

398,100 

478,100 

67,400 

95,200 

63,100 

Failure  Criterion 
(x  Initial  LVDT) 

Load 

Frequency 

(Hz) 

cn 

m 

m 

ro 

Max.  Load 
Level  (lbs) 

12,910 

12,929 

12,753 

18,971 

18,428 

18,268 

12,751 

12,570 

12,742 

18,434 

18,676 

18,169 

Max-Min 
Stress  (ksi) 

±5.426 

±5.426 

±5.426 

±7.597 

±7.597 

±7.597 

±5.426 

±5.426 

±5.426 

±7.597 

±7.597 

±7.597 

RMS 

(ksi) 

1.25 

1.25 

1.25 

1.75 

1.75 

1 

1.25 

1.25 

1.75 

1.75 

1.75 

Identification 

No. 

7-0750-10 

7-0750-13 

1  . . 

8-0750-7 

1 _ 

7-0750-11 

7-0750-14 

8-0750-5 

r- 

1 

o 

o 

1 

VO 

6-0750-8 

7-0750-4 

6-0750-6 

7-0750-1 

7-0750-3 

Lay¬ 

up 

O 

5 

0/90° 
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TABLE  6  -  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  FRP  JOINTS 


Q 

73 

e 

o 

•  M 

-3 
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< 
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S 

u 
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1 

1 

1 

1 

j 

! 

1 

1 

1 
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1 

1 

1 

1 

1 

I 

1 

H 

g 

B 
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1 

s 

1 

1 

B 

s 

s 

s 

B 

s 

1 

1 

1 

3 
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u 

o 

o 
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(N 

o 

o 

a\ 
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QO 
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o 

in 

O 

o 

in 

m 

CN 

o 

o 

m 

o' 

o 

o 

o 
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o 

o 

CN 
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00 

o 
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o 

CN 

CN 

o 

00 

ON 

o 

VO 
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(x  Initial  LVDT) 

tT 

Load 
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(Hz) 

CN 

- 

- 

- 

1 

1 

- 

m 

cn 

Max.  Load 
Level  (lbs) 

in 

CN| 

trT 

xn 

m 

in 

in 

cn 

r-" 

ON 

in 

m 

i> 

r-" 

cn 

cP 

o 

r- 

o" 

VO 

O 

00^ 

CN 

VO 

00 

CN 

r- 

m 

cc 

1-H 

rr 

VO 
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CN 

CN 

cn 

cn 

in 

in 

in 

in 

c 

o 
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-o 

NH 

O 

1 

o 

CN 

1 

c 
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1 

00 

1 

< 
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00 

1 

< 
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T3 

(U 

-o 
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TABLE  6  (con’t)  -  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  FRP  JOINTS 


Additional  Data 

Temp. 

I 

f 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

LVDT 

t 

1 

S 

s 

s 

s 

S 
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1 

1 

1 

s 

s 

s 

S 

Cycles  to 
Failure 

589,400 

20,000,000  + 

23,200 

899,100 

12,600 

8,720 

1,810 

2,530 

o 

o^ 

00 

O 

o 

00 

305 

320 

Failure  Criterion 
(x  Initial  LVDT) 

rr 

cn 

Load 

Frequency 

(Hz) 

- 

- 

T-H 

- 

- 

- 

- 

- 

- 

CO 

0.5 

Max.  Load 
Level  (lbs) 

7,607 

7,019 

7,307 

7,319 

8,617 

8,409 

9,623 

9,464 

11,795 

11,695 

17,774 

17,793 

Max.  Stress 
Level  (ksi) 

m 

m 

m 

ro 

3.5 

3.5 

iT) 

7.5 

7.5 

Identification 

Name 

JB-7 

oo 

1 

OQ 

1-^ 

JB-13 

JB-14 

JB-12 

JB-9 

JB-  10 

JB-5 

JB-6 

JB-3 

JB-4 

Joint 

Type 

Bolted 
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Indicates  test  was  suspended  without  failure. 


TABLE  6  (con’t)  -  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  FRP  JOINTS 


S 
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TABLE  7  -  VARIABLE  AMPLITUDE  FATIGUE  TESTS  OF  FRP  JOINTS 
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-  17  2.25  ±9.765  24,664  3  4  67,400 


TABLE  8  -  COMPARISON  OF  RAYLEIGH  FATIGUE  LIVES  WITH  TEST  RESULTS 
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